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This is a short review of the theoretical work on the two-dimensional Hubbard model performed 
in Sherbrooke in the last few years. It is written on the occasion of the twentieth anniversary of 
the discovery of high-temperature superconductivity. We discuss several approaches, how they were 
benchmarked and how they agree sufficiently with each other that we can trust that the results are 
accurate solutions of the Hubbard model. Then comparisons are made with experiment. We show 
that the Hubbard model does exhibit d-wave superconductivity and antiferromagnetism essentially 
where they are observed for both hole and electron-doped cuprates. We also show that the pseudo- 
gap phenomenon comes out of these calculations. In the case of electron-doped high temperature 
superconductors, comparisons with angle-resolved photoemission experiments are nearly quantita- 
tive. The value of the pseudogap temperature observed for these compounds in recent photoemission 
experiments has been predicted by theory before it was observed experimentally. Additional exper- 
imental confirmation would be useful. The theoretical methods that are surveyed include mostly 
the Two-Particle Self-Consistent Approach, Variational Cluster Perturbation Theory (or variational 
cluster approximation), and Cellular Dynamical Mean-Field Theory. 



Contents 



I. Introduction 1 

II. Methodology 2 

A. Weak coupling approach 3 

1. Two-Particle Self-consistent approach 
(TPSC) ' 3 

2. Benchmarks for TPSC 5 

B. Strong-coupling approaches: Quantum 
clusters 7 

1. Cluster perturbation theory 8 

2. Self-energy functional approach 8 

3. Variational cluster perturbation theory, or 
variational cluster approximation 9 

4. Cellular dynamical mean-field theory 10 

5. The Dynamical cluster approximation 10 

6. Benchmarks for quantum cluster 
approaches 1 1 

III. Results and concordance between different 
methods 12 

A. Weak and strong-coupling pseudogap 12 

1. CPT 13 

2. CDMFT and DCA 14 

3. TPSC (including analytical results) 15 

4. Weak- and strong-coupling pseudogaps 17 

B. d-wave superconductivity 18 

1. Zero-temperature phase diagram 18 

2. Instability of the normal phase 20 

IV. Comparisons with experiment 21 

A. ARPES spectrum, an overview 21 

B. The pseudogap in electron-doped cuprates 22 

C. The phase diagram for high-temperature 
superconductors 25 



1. Competition between antiferromagnetism 



and superconductivity 25 
2. Anomalous superconductivity near the 

Mott transition 27 

V. Conclusion, open problems 27 

Acknowledgments 29 

A. List of acronyms 29 

References 30 



I. INTRODUCTION 

In the first days of the discovery of high-temperature 
superconductivity, Anderson 1 suggested that the two- 
dimensional Hubbard model held the key to the phe- 
nomenon. Despite its apparent simplicity, the two- 
dimensional Hubbard model is a formidable challenge for 
theorists. The dimension is not low enough that an exact 
solution is available, as in one dimension. The dimension 
is not high enough that some mean-field theory, like Dy- 
namical Mean Field Theory 2,3 (DMFT), valid in infinite 
dimension, can come to the rescue. In two dimensions, 
both quantum and thermal fluctuations are important. 
In addition, as we shall see, it turns out that the real 
materials are in a situation where both potential and 
kinetic energy are comparable. We cannot begin with 
the wave picture (kinetic energy dominated, or so-called 
"weak coupling") and do perturbation theory, and we 
cannot begin from the particle picture (potential energy 
dominated, or so-called "strong coupling" ) and do per- 
turbation theory. In fact, even if one starts from the 
wave picture, perturbation theory is not trivial in two 
dimensions, as we shall see. Variational approaches on 
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FIG. 1: The d x i_ y i and extended s-wave susceptibilities ob- 
tained from QMC simulations for U = 4t and a 4 x 4 lattice. 
The solid lines are the non-interacting results. From Ref. 5. 
The low temperature downturn of d-wave however seems to 
come from a mistreatment of the sign problem (D.J. Scalapino 
private communication) . 



the ground state have been proposed, 4 but even if they 
capture key aspects of the ground state, they say little 
about one-particle excitations. 

Even before the discovery of high-temperature super- 
conductivity, it was suggested that antiferromagnetic 
fluctuations present in the Hubbard model could lead to 
d-wave superconductivity, 6-8 a sort of generalization of 
the Kohn-Luttinger mechanism 9 analogous to the super- 
fluidity mediated by ferromagnetic spin fluctuations in 
3 He. 10 Nevertheless, early Quantum Monte Carlo (QMC) 
simulations 5 gave rather discouraging results, as illus- 
trated in Fig. 1. In QMC, low temperatures are in- 
accessible because of the sign problem. At accessible 
temperatures, the d-wave pair susceptibility is smaller 
than the non-interacting one, instead of diverging. Since 
the observed phenomenon appears at temperatures that 
are about ten times smaller than what is accessible with 
QMC, the problem was left open. Detailed analysis of 
the irreducible vertex 45 deduced from QMC did suggest 
the importance of d-wave pairing, but other numerical 
work 12 concluded that long-range d-wave order is ab- 
sent, despite the fact that slave-boson approaches 13 ' 14 
and many subsequent work suggested otherwise. The 
situation on the numerical side is changing since more 
recent variational, 4 Dynamical Cluster Approximation 15 
and exact diagonalization 16 results now point towards the 
existence of d-wave superconductivity in the Hubbard 
model. Even more recently, new numerical approaches 
are making an even more convincing case. 17-19 

After twenty years, we should be as quantitative as 
possible. How should we proceed to investigate a model 
without a small parameter? We will try to follow this 
path: (1) Identify important physical principles and laws 
to constrain non-perturbative approximation schemes, 
starting from both weak (kinetic energy dominated) 
and strong (potential energy dominated) coupling. (2) 
Benchmark the various approaches as much as possi- 
ble against exact (or numerically accurate) results. (3) 
Check that weak and strong coupling approaches agree 
at intermediate coupling. (4) Compare with experiment. 



In brief, we are trying to answer the question, "Is 
the Hubbard model rich enough to contain the essential 
physics of the cuprates, both hole and electron doped?" 
The answer is made possible by new theoretical ap- 
proaches, increased computing power, and the reassur- 
ance that theoretical approaches, numerical and analyti- 
cal, give consistent results at intermediate coupling even 
if the starting points are very different. 

This paper is a review of the work we have done in 
Sherbrooke on this subject. In the short space provided, 
this review will not cover all of our work. Needless to 
say, we will be unfair to the work of many other groups, 
even though we will try to refer to the work of others 
that is directly relevant to ours. We do not wish to make 
priority claims and we apologize to the authors that may 
feel unfairly treated. 

Section II will introduce the methodology: First a 
method that is valid at weak to intermediate coupling, 
the Two-Particle Self- Consistent approach (TPSC), and 
then various quantum cluster methods that are better 
at strong coupling, namely Cluster Perturbation Theory 
(CPT), the Variational Cluster Approximation (VCA) 
also known as Variational Cluster Perturbation Theory 
(VCPT), and Cellular Dynamical Mean Field Theory 
(CDMFT) with a brief mention of the Dynamical Cluster 
Approximation (DC A). In all cases, we will mention the 
main comparisons with exact or numerically accurate re- 
sults that have been used to benchmark the approaches. 
In Sect. Ill we give some of the results, mostly on the 
pscudogap and the phase diagram of high-temperature 
superconductors. More importantly perhaps, we show 
the consistency of the results obtained by both weak- and 
strong-coupling approaches when they are used at inter- 
mediate coupling. Finally, we compare with experiment 
in section IV. 



II. METHODOLOGY 

We consider the Hubbard model 

i,j,cr i 

where c\ a (c i(T ) are creation and annihilation operators 
for electrons of spin a, n ia — c\ a c ia is the density of 
spin a electrons, tij — £!•» is the hopping amplitude, 
and U is the on-site Coulomb repulsion. In general, we 
write t, t' , t" respectively for the first-, second- and third- 
nearest neighbor hopping amplitudes. 

In the following subsections, we first discuss how to ap- 
proach the problem from the weak coupling perspective 
and then from the strong coupling point of view. The ap- 
proaches that we will use in the end are non-perturbative, 
but in general they are more accurate either at weak or 
strong coupling. 
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A. Weak coupling approach 

Even at weak coupling, the Hubbard model presents 
difficulties specific to two dimensions. The time-honored 
Random Phase Approximation (RPA) has the advantage 
of satisfying conservation laws, but it violates the Pauli 
principle and the Mermin-Wagncr-Hohcnberg-Coleman 
(or Mermin- Wagner, for short) theorem. This theorem 
states that a continuous symmetry cannot be broken 
at finite temperature in two dimensions. RPA gives a 
finite-temperature phase transition. The Pauli principle 
means, in particular, that (n^n^) = (n^) in a model 
with only one orbital per site. This is violated by RPA 
since it can be satisfied only if all possible exchanges 
of electron lines are allowed (more on this in the fol- 
lowing section). Since the square of the density at a 
given site is given by ((n^ + n i \ l ) 2 ) = 2{n-\n-\) + 2{m i n{) 1 
violating the Pauli condition (n^n^) = (n^) will in 
general lead to large errors in double occupancy, a 
key quantity in the Hubbard model since it is propor- 
tional to the potential energy. Another popular ap- 
proach is the Moriya 20 self-consistent spin-fluctuation 
approach 21 that uses a Hubbard-Stratonovich transfor- 
mation and a {<p 4 ) ~ <p 2 {<p 2 ) factorization. This satis- 
fies the Mermin- Wagner theorem but, unfortunately, vi- 
olates the Pauli principle and introduces an unknown 
mode-coupling constant as well as an unknown renor- 
malized U in the second-order term. The conserving 
approximation known as Fluctuation Exchange (FLEX) 
approximation 22 is an Eliashbcrg-type theory that is con- 
serving but violates the Pauli principle, assumes a Migdal 
theorem and does not reproduce the pseudogap phe- 
nomenon observed in QMC. More detailed criticism of 
this and other approaches may be found in Refs. 23,24. 
Finally, the renormalization group 25-29 has the great ad- 
vantage of being an unbiased method to look for instabil- 
ities towards various ordered phases. However, it is quite 
difficult to implement in two dimensions because of the 
proliferation of coupling constants, and, to our knowl- 
edge, no one has yet implemented a two-loop calcula- 
tion without introducing additional approximations. 30 ' 31 
Such a two-loop calculation is necessary to observe the 
pseudogap phenomenon. 



1. Two-Particle Self- consistent approach (TPSC) 

The TPSC approach, originally proposed by Vilk, 
Tremblay and collaborators, 33 ' 37 aims at capturing non- 
perturbative effects. It does not use perturbation theory 
or, if you want, it drops diagrammatic expansions. In- 
stead, it is based on imposing constraints and sum rules: 
the theory should satisfy (a) the spin and charge conser- 
vation laws (b) the Pauli principle in the form {n^n^) = 
(riii) (c) the local- moment and the local-density sum 
rules. Without any further explicit constraint, we find 
that the theory satisfies the Mermin- Wagner theorem, 
that it satisfies consistency between one- and two-particle 



quantities in the sense that ±Tr(£G) = U{n^n{) and fi- 
nally that the theory contains the physics of Kanamori- 
Briickner screening (in other words, scattering between 
electrons and holes includes T-matrix quantum fluctua- 
tion effects beyond the Born approximation). 

Several derivations of our approach have been 
given, 32,33 including a quite formal one 34 based on the 
functional derivative Baym-Kadanoff approach. 35 Here 
we only give an outline 36 of the approach with a more 
phenomenological outlook. We proceed in two steps. In 
the first step (in our earlier work sometimes called zeroth 
step), the self-energy is obtained by a Hartree-Fock-type 
factorization of the four-point function with the addi- 
tional constraint that the factorization is exact when all 
space-time coordinates coincide. 157 Functional differenti- 
ation, as in the Baym-Kadanoff approach 35 , then leads 
to a momentum- and frequency-independent irreducible 
particle-hole vertex for the spin channel that satisfies 37 
U sp = U{n^n\) /{(n-\)(ni)). The local moment sum rule 
and the Pauli principle in the form (n 2 ) = (n CT ) then 
determine double occupancy and U sp . The irreducible 
vertex for the charge channel is too complicated to be 
computed exactly, so it is assumed to be constant and 
its value is found from the Pauli principle and the local 
charge fluctuation sum rule. To be more specific, let us 
use the notation, q = (q,iq n ) and k — (k,ik n ) with iq n 
and ik n respectively bosonic and fermionic Matsubara 
frequencies. We work in units where ks, ft, and lattice 
spacing are all unity. The spin and charge susceptibilities 
now take the form 



and 



x C h(q) = xo(qY 



-u ch 



(2) 



(3) 



with xo computed with the Green function that con- 
tains the self-energy whose functional differentiation gave 
the vertices. This self-energy is constant, corresponding 
to the Hartree-Fock-type factorization. 158 The suscepti- 
bilities thus satisfy conservation laws 35 . One enforces the 
Pauli principle (n 2 ) = (n a ) implicit in the following two 
sum rules, 



T 

N 

T 



Xsp(q) = ((«t - n if) = n - 2 (™T n l) 



(4) 



Xch(q) = ((n T + ra x ) 2 ) - n 2 = n + 2(n T n x ) - n 2 



where n is the density. The above equations, in addition 



to' 



37 



U s 



C/(n T n x ) 



(«T>(-l)' 

suffice to determine the constant vertices U sp and U c h- 



(5) 
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Once the two-particle quantities have been found as 
above, the next step of the approach of Ref. 23, consists in 
improving the approximation for the single-particle self- 
energy by starting from an exact expression where the 
high-frequency Hartree-Fock behavior is explicitly fac- 
tored out. One then substitutes in the exact expression 
the irreducible low- frequency vertices U sp and U c h as well 

as G^\k + q) and Xsp(q),Xch(q) computed above. The 
exact form for the self-energy expression can however be 
obtained either in the longitudinal or in the transverse 
channel. To satisfy crossing symmetry of the fully re- 
ducible vertex appearing in the general expression and to 
preserve consistency between one- and two-particle quan- 
tities, one averages the two possibilities to obtain 36 

+ | ^ E ^ U spXs P (q) + U chXch (q)} G^(k + q). (6) 
Q 

The resulting self-energy (fc) on the left hand- 
side is at the next level of approximation so it dif- 
fers from the self-energy entering the right-hand side. 
One can verify that the longitudinal spin fluctuations 
contribute an amount U{n-\n{)/A to the consistency 
condition 24 iTr(£( 2 )GW) = U{n^ ni ) and that each of 
the two transverse spin components as well as the charge 
fluctuations also each contribute U{n^n\ l )/A. In addi- 
tion, one verifies numerically that the exact sum rule 23 
— J du' lm[T, a (k,uj')]/'K = C/ 2 n_ cr (l — n_ CT ) determining 
the high-frequency behavior is satisfied to a high degree 
of accuracy. 

The theory also has a consistency check. Indeed, the 
exact expression for consistency between one- and two- 
particle quantities should be written with G^ 2 -* given by 



(G*- 1 )^) = (G- 1 )^) -£< 2 ) instead of with G {1 \ In other 
words iTr(S( 2 )G( 2 )) = U^n^ should be satisfied in- 
stead of iTr(£( 2 )G (1) ) = U(n^n{), which is exactly satis- 
fied here. We find through QMC benchmarks that when 
the left- and right-hand side of the last equation differ 
only by a few percent, then the theory is accurate. 

To obtain the thermodynamics, one finds the entropy 
by integrating 1/T times the specific heat (dE/dT) so 
that we know F = E - TS. There are other ways to ob- 
tain the thermodynamics and one looks for consistency 
between these. 38 We will not discuss thermodynamic as- 
pects in the present review. 

At weak coupling in the repulsive model the particle- 
hole channel is the one that is influenced directly 
Correlations in crossed channels, such as pairing sus- 
ceptibilities, are induced indirectly and are harder to 
evaluate. This simply reflects the fact the simplest 
Hartree-Fock factorization of the Hubbard model does 
not lead to a d-wave order parameter (even though 
Hartree-Fock factorization of its strong-coupling version 
does). The d a .2_ y 2-wave susceptibility is defined by 

Xd = Jq dr(r r A(t)A^) with the d-wave order param- 
eter equal to A* = J2i S 7 3(7) c !t c !+71 tnc sum over 7 
being over nearest-neighbors, with 3(7) = ±1/2 depend- 
ing on whether 7 is a neighbor along the x or the y axis. 
Briefly speaking, 39 ' 40 to extend TPSC to compute pair- 
ing susceptibility, we begin from the Schwinger-Martin- 
Kadanoff-Baym formalism with both diagonal 23 ' 34 and 
off-diagonal 41 source fields. The self-energy is expressed 
in terms of spin and charge fluctuations and the irre- 
ducible vertex entering the Bcthe-Salpeter equation for 
the pairing susceptibility is obtained from functional dif- 
ferentiation. The final expression for the d-wave suscep- 
tibility is, 



xM = o,iq n = 0) = ^Yl (g 2 d ^)Gf\-k)Gf\k)) u - 5>(k)G; 2) (-fc)Gf } (fc) 

k ^ ' k,k' 

x f 1 T TT^~lh' — h\ + VTTTr 1 <v G?(-k>)GV(k')g d (k!). (7) 

\1 - 2 U s P Xo{k - k) 1 + zUchXoik '- k)J ' 



In the above expression, gd(k) is the form factor for the 
gap symmetry, while k and k! stand for both wave-vector 
and fermionic Matsubara frequencies on a square-lattice 
with N sites at temperature T. The spin and charge sus- 
ceptibilities take the form x7p(l) — Xoiq)' 1 — \U sp and 
Xch(q) = Xo(q)^ 1 + \U c h with xo computed with the 
Green function G^ that contains the self-energy whose 
functional differentiation gave the spin and charge ver- 
tices. The values of U sp , U c h and (n-j-nj.) are obtained as 
described above. In the pseudogap regime, one cannot 
use U sp — U{n^n\)/{{n^){n\)). Instead, 23 one uses the 



I 

local-moment sum rule with the zero temperature value 
of (n^ni) obtained by the method of Ref. 42 that agrees 

very well with QMC calculations at all values of U. Also, 

(2) 

GJr contains self-energy effects coming from spin and 
charge fluctuations, as described above. 34 ' 36 

The same principles and methodology can be applied 
for the attractive Hubbard model. 39 ' 41,43 In that case, the 
dominant channel is the s-wave pairing channel. Cor- 
relations in the crossed channel, namely the spin and 
charge susceptibilities, can also be obtained mutatis mu- 
tandi along the lines of the previous paragraph. 
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FIG. 2: Comparisons between the QMC simulations (sym- 
bols) and TPSC (solid lines) for the filling dependence of the 
double occupancy. The results are for T — t/6 as a function of 
filling and for various values of U expect for U = At where the 
dashed line shows the results of our theory at the crossover 
temperature T = T X . From Ref. 23. 



2. Benchmarks for TPSC 

To test any non-perturbative approach, we need reli- 
able benchmarks. Quantum Monte Carlo (QMC) sim- 
ulations provide such benchmarks. The results of such 
numerical calculations are unbiased and they can be ob- 
tained on much larger system sizes than any other simula- 
tion method. The statistical uncertainty can be made as 
small as required. The drawback of QMC is that the sign 
problem renders calculations impossible at temperatures 
low enough to reach those that are relevant for d-wave 
superconductivity. Nevertheless, QMC can be performed 
in regimes that are non-trivial enough to allow us to elim- 
inate some theories on the grounds that they give quali- 
tatively incorrect results. Comparisons with QMC allow 
us to estimate the accuracy of the theory. An approach 
like TPSC can then be extended to regimes where QMC 
is unavailable with the confidence provided by agreement 
between both approaches in regimes where both can be 
performed. 

In order to be concise, details are left to figure cap- 
tions. Let us first focus on quantities related to spin and 
charge fluctuations. The symbols on the figures refer to 
QMC results while the solid lines come from TPSC cal- 
culations. Fig. 2 shows double occupancy, a quantity 
that plays a very important role in the Hubbard model 
in general and in TPSC in particular. That quantity is 
shown as a function of filling for various values of U at 
inverse temperature (3 = 6. Fig. 3 displays the spin and 
charge structure factors in a regime where size effects are 
not important. Clearly the results are non-perturbative 
given the large difference between the spin and charge 
structure factors, which are plotted here in units where 
they are equal at U = 0. In Fig. 4 we exhibit the static 
structure factor at half-filling as a function of tempera- 
ture. Below the crossover temperature Tx, there is an 
important size dependence in the QMC results. The 
TPSC calculation, represented by a solid line, is done 
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FIG. 3: Wave vector (q) dependence of the spin and charge 
structure factors for different sets of parameters. Solid lines 
are from TPSC and symbols are our QMC data. Monte Carlo 
data for n = 1 and U = 8t are for 6x6 clusters and T = 0.5t; 
all other data are for 8x8 clusters and T = 0.2t. Error bars 
are shown only when significant. From Ref. 37. 




FIG. 4: Temperature dependence of S ap (iT,n) at half-filling 
n = 1. The solid line is from TPSC and symbols are Monte 
Carlo data from Ref. 44. Taken from Ref. 37. 



for an infinite system. We see that the mean-field finite 
transition temperature Tmf is replaced by a crossover 
temperature Tx at which the correlations enter an ex- 
ponential growth regime. One can show analytically 23,37 
that the correlation length becomes infinite only at zero 
temperature, thus satisfying the Mermin- Wagner theo- 
rem. The QMC results approach the TPSC results as the 
system size grows. Nevertheless, TPSC is in the N = oo 
universality class 46 contrary to the Hubbard model for 
which N = 3, so one expects quantitative differences to 
increase as the correlation length becomes larger. It is 
important to note that Tx does not coincide with the 
mean- field transition temperature Tmf- This is because 
of Kanamori-Brueckner screening 37,47 that manifests it- 
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FIG. 5: Comparisons between Monte Carlo simulations 
(BW), FLEX calculations and TPSC for the spin suscepti- 
bility at Q = (tv, 7t) as a function of temperature at zero Mat- 
subara frequency. The filled circles (BWS) are from Ref. 45. 
Taken from Ref. 23. 
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FIG. 6: Single particle spectral weight A(k,cj) for U = At, 
j3 = 5/t, n = 1, and all independent wave vectors k of an 
8x8 lattice. Results obtained from Maximum Entropy in- 
version of QMC data on the left panel and many-body TPSC 
calculations with Eq.( 6) on the middle panel and with FLEX 
on the right panel. From Ref. 36. 



self in the difference between U sp and the bare U . 



Be- 



low Tx, the main contribution to the static spin struc- 
ture factor in Fig. 4 comes from the zero-Matsubara fre- 
quency component of the spin susceptibility. This is the 
so-called renormalized classical regime where the charac- 
teristic spin fluctuation frequency u sp is much less than 
temperature. Even at temperatures higher than that, 
TPSC agrees with QMC calculation much better than 
other methods, as shown in Fig. 5. 

Below the crossover temperature to the renormalized 
classical regime, a pseudogap develops in the single- 
particle spectral weight. This is illustrated in Fig. 6. 36 
Eliashberg-type approaches such as FLEX do not show 
the pseudogap present in QMC. The size dependence of 
the results is also quite close in TPSC and in QMC, as 
shown in Fig. 7. 

The d-wave susceptibility 40 shown in Fig. 8 again 
clearly demonstrates the agreement between TPSC and 
QMC. In particular, the dome shape dependence of the 
QMC results is reproduced to within a few percent. We 
will see in Sec. Ill how one understands the dome shape 
and the fact that the d-wave susceptibility of the inter- 
acting system is smaller than that of the non-interacting 




FIG. 7: Size dependent results for various types of calcula- 
tions for U = At, (3 = 5/t, n = 1, k = (0,7r), L = 4,6,8, 10. 
Upper panels show A(k, ui) extracted from Maximum Entropy 
on G(t) shown on the corresponding lower panels, (a) QMC. 
(b) TPSC using Eq. (6). (c) FLEX. From Ref. 36. 
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FIG. 8: Comparisons between the d x 2_ y 2 susceptibility ob- 
tained from QMC simulations (symbols) and from the TPSC 
approach (lines) in the two-dimensional Hubbard model. 
Both calculations are for U — At, a 6 x 6 lattice. QMC er- 
ror bars are smaller than the symbols. Analytical results are 
joined by solid lines. The size dependence of the results is 
small at these temperatures. The U = case is also shown 
at P = A/t as the upper line. The inset compares QMC and 
FLEX at U = A,P = A/t. From Ref. 40. 



one in this temperature range. 

To conclude this section, we quickly mention a few 
other results obtained with TPSC. Fig. 9 contrasts the 
crossover phase diagram obtained for the Hubbard model 
at the van Hove filling 48 with the results of a renormaliza- 
tion group calculation. 28 The difference occurring in the 
ferromagnetic region is discussed in detail in Ref. 48. Fi- 
nally, we point out various comparisons for the attractive 
Hubbard model. Fig. 10 shows the static s-wave pairing 
susceptibility, Fig. 11 the chemical potential and the oc- 
cupation number, and finally Fig. 12 the local density of 
states and the single-particle spectral weight at a given 
wave vector. 



FIG. 9: The crossover diagram as a function of next-nearest- 
neighbor hopping t' from TPSC (left) and from a temperature 
cutoff renormalization group technique from Ref. 28 (right). 
The corresponding Van Hove filling is indicated on the up- 
per horizontal axis. Crossover lines for magnetic instabilities 
near the antiferromagnetic and ferromagnetic wave vectors 
are represented by filled symbols while open symbols indicate 
instability towards d x 2_ y 2-wave superconducting. The solid 
and dashed lines below the empty symbols show, respectively 
for U — 3t and U = 6t, where the antiferromagnetic crossover 
temperature would have been in the absence of the super- 
conducting instability. The largest system size used for this 
calculation is 2048 x 2048. From Ref. 48. 
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FIG. 10: TPSC s-wave paring structure factor S(q, r = 0) 
(filled triangles) and QMC 5*(q, r = 0) (open circles) for U = 
—At and various temperatures (a) at n — 0.5 and (b) at n = 
0.8 on a 8 x 8 lattice. The dashed lines are to guide the eye. 
From Ref. 43. 



B. Strong-coupling approaches: Quantum clusters 

DMFT 3,51 has been extremely successful in helping us 
understand the Mott transition, the key physical phe- 
nomenon that manifests itself at strong coupling. How- 
ever, in high dimension where this theory becomes ex- 
act, spatial fluctuations associated with incipient order 
do not manifest themselves in the self-energy. In low di- 
mension, this is not the case. The self-energy has strong 
momentum dependence, as clearly shown experimentally 
in the high-temperature superconductors, and theoreti- 
cally in the TPSC approach, a subject we shall discuss 
again below. It is thus necessary to go beyond DMFT 
by studying clusters instead of a single Anderson impu- 



FIG. 11: Left: chemical potential shifts — fio (open dia- 
monds) and /i' 2 ' — fiQ (open squares) with the results of QMC 
calculations (open circles) for U = —At. Right: The momen- 
tum dependent occupation number n(k). Circles: QMC cal- 
culations from Ref. 49. The solid curve: TPSC. The dashed 
curve obtained by replacing U pp by U in the self-energy with 
all the rest unchanged. The long-dash line is the result of a 
self-consistent T-matrix calculation, and the dot-dash line the 
result of second-order perturbation theory. From Ref. 43. 
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FIG. 12: Comparisons of local density of states and single- 
particle spectral weight from TPSC (solid lines) and QMC 
(dashed lines) on a 8 x 8 lattice. QMC data for the density 
of states taken from Ref. 50. Figures from Ref. 43. 



rity as done in DMFT. The simplest cluster approach 
that includes strong-coupling effects and momentum de- 
pendence is Cluster Perturbation Theory (CPT). 52,53 In 
this approach, an infinite set of disconnected clusters are 
solved exactly and then connected to each other using 
strong-coupling perturbation theory. Although the re- 
sulting theory turns out to give the exact result in the 
[7 = case, its derivation clearly shows that one expects 
reliable results mostly at strong coupling. This approach 
does not include the self-consistent effects contained in 
DMFT. Self-consistency or clusters was suggested in Ref. 
2,54 and a causal approach was first implemented within 
DCA, 55 where a momentum-space cluster is connected 
to a self-consistent momentum-space medium. In our 
opinion, the best framework to understand all other clus- 
ter methods is the Self-Energy Functional approach of 
Potthoff. 56,57 The form of the lattice Green function ob- 
tained in this approach is the same as that obtained in 
CPT, clearly exhibiting that such an approach is better 
at strong-coupling, even though results often extrapolate 
correctly to weak coupling. Amongst the special cases of 
this approach, the Variational Cluster Approach (VCA), 
or Variational Cluster Perturbation Theory (VCPT) 57 is 
the one closest to the original approach. In a variant, 
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Cellular Dynamical Mean Field Theory 58 (CDMFT), a 
cluster is embedded in a self-consistent medium instead 
of a single Anderson impurity as in DMFT (even though 
the latter approach is accurate in many realistic cases, 
especially in three dimensions). The strong-coupling as- 
pects of CDMFT come out clearly in Refs. 59,60. A de- 
tailed review of quantum cluster methods has appeared 
in Ref. 61. 



1. Cluster perturbation theory 

Even though CPT does not have the self-consistency 
present in DMFT type approaches, at fixed computing 
resources it allows for the best momentum resolution. 
This is particularly important for the ARPES pseudogap 
in electron-doped cuprates that has quite a detailed mo- 
mentum space structure, and for d-wave superconducting 
correlations where the zero temperature pair correlation 
length may extend well beyond near-neighbor sites. CPT 
was developed by Gros 52 and Senechal 53 independently. 
This approach can be viewed as the first term of a sys- 
tematic expansion around strong coupling. 62 Let us write 
the hopping matrix elements in the form 



f mn _ .(c) r 



_i_ umn 



(8) 



where m and n label the different clusters, and /i, v label 

(c) 

the sites within a cluster. Then t^l labels all the hopping 
matrix elements within a cluster and the above equation 
defines V™. 

We pause to introduce the notation that will be used 
throughout for quantum cluster methods. We follow the 
review article Ref. 61. In reciprocal space, any wave vec- 
tor k in the Brillouin zone may be written as k =k + K 
where both k and k are continuous in the infinite size 
limit, except that k is defined only in the reduced Bril- 
louin zone that corresponds to the superlattice. On the 
other hand, K is discrete and denotes reciprocal lattice 
vectors of the superlattice. By analogy, any position r 
in position space can be written as f + R where R is for 
positions within clusters while f labels the origins of the 
clusters, an infinite number of them. Hence, Fourier's 
theorem allows one to define functions of k, k or K 
that contain the same information as functions of, re- 
spectively, r, f or R. Also, we have K f = 2im where n 
is an integer. Sites within a cluster are labelled by greek 
letters so that the position of site \i within a cluster is 
R M , while clusters are labelled by Latin letters so that 
the origin of cluster m is at f m . 

Returning to CPT, the Green function for the whole 
system is given by 



6 -1 (k,«)l =\&v>- 1 (z)-v$) 



(9) 



where hats denote matrices in cluster site indices and z 
is the complex frequency. At this level of approximation, 
the CPT Green function has the same structure as in 



the Hubbard I approximation except that it pertains to 
a cluster instead of a single site. Since G^~ x {z) = z + 
jj-ito-EW and G(°)- 1 (k,2)=z + /i-iW- V(k), the 
Green function (9) may also be written as 



G- 1 (lL,z) = G^- 1 (^,z)-t^(z 



(10) 



This form allows a different physical interpretation of the 
approach. In the above expression, the self-energy of the 
lattice is approximated by the self-energy of the cluster. 
The latter in real space spans only the size of the cluster. 

We still need an expression to extend the above result 
to the lattice in a translationally invariant way. This is 
done by defining the following residual Fourier transform: 



1 N c 

G CPT (k,z) = — £e ik '( R -- R ">G^(M. 



(11) 



Notice that G M „(k, z) may be replaced by G M „(k, z) in 
the above equation since V^k + K) = V"(k). 

2. Self-energy functional approach 

The self-energy functional approach, devised by 
Potthoff 57 allows one to consider various cluster schemes 
from a unified point of view. It begins with f2t[G], a 
functional of the Green function 

Q t [G] = $[G] - Tr((Gj£ - G- X )G) + Trln(-G). (12) 

The Luttinger Ward functional $[G] entering this equa- 
tion is the sum of connected vacuum skeleton diagrams. 
A diagram-free definition of this functional is also given 
in Ref. 63. For our purposes, what is important is that 
(1) The functional derivative of $[G] is the self-energy 



6$[G] 
5G 



= £ 



(13) 



and (2) it is a universal functional of G in the following 
sense: whatever the form of the one-body Hamiltonian, 
it depends only on the interaction and, functionnally, it 
has the same dependence on G. The dependence of the 
functional Q t [G] on the one-body part of the Hamiltonian 
is denoted by the subscript t and it comes only through 
Gqj. 1 appearing on the right-hand side of Eq. (12). 

The functional Qt [G] has the important property that 
it is stationary when G takes the value prescribed by 
Dyson's equation. Indeed, given the last two equations, 
the Euler equation takes the form 



sn t [G] 

5G 



S - Go/ + G- 1 = 0. 



(14) 



This is a dynamic variational principle since it involves 
the frequency appearing in the Green function, in other 
words excited states are involved in the variation. At this 
stationary point, and only there, fit[G] is equal to the 



9 



' »o ua ma 
p»1p »k> »K 

i • o • a • 01 
p »lp •to *\ 

> • oi • a • a 
g | p • g » 



o • o • 
• o • o 


o • o • 
• o • o 


o • o • 
• o • o 


• • 

• o • o 


• o • 
• o • o 


• • 

• o • o 


o • o • 
• • o 


o • o • 

• 0*0 


o • o • 
• o • o 



o 
o o 


o 
o o 


o 
o o 


o 
o o 


o 
o o 


o 
o o 


o 
o o 


o 
o o 


o 
o o 




FIG. 13: Various tilings used in quantum cluster approaches. 
In these examples the grey and white sites are inequivalent 
since an antiferromagnetic order is possible. 



grand potential. Contrary to Ritz's variational principle, 
this last equation does not tell us whether f2t[G] is a 
minimum or a maximum or a saddle point there. 

There are various ways to use the stationarity prop- 
erty that we described above. The most common one, 
is to approximate <&[G] by a finite set of diagrams. 
This is how one obtains the Hartree-Fock, the FLEX 
approximation 22 or other so-called thermodynamically 
consistent theories. This is what Potthoff calls a type II 
approximation strategy. 64 A type I approximation sim- 
plifies the Euler equation itself. In a type III approxi- 
mation, one uses the exact form of $[G] but only on a 
limited domain of trial Green functions. 

Following Potthoff, we adopt the type III approxima- 
tion on a functional of the self-energy instead of on a 
functional of the Green function. Suppose we can locally 
invert Eq. (13) for the self-energy to write G as a func- 
tional of E. We can use this result to write, 



fit[E] = F[E] - Trln(- 



£). 



where we defined 



F[E] = *[G] -Tr(EG). 



(15) 



(16) 



and where it is implicit that G = G[S] is now a func- 
tional of E. F[E], along with the expression (13) for the 
derivative of the Luttinger-Ward functional, define the 
Legendre transform of the Luttinger-Ward functional. It 
is easy to verify that 



6F[E] c5$[G]<5G[E] £G[E] 



<5E 



SG 5E 



5£ 



G 



-G (17) 



hence, fit [E] is stationary with respect to E when Dyson's 
equation is satisfied 



Mjt[gj 
5E 



-G+(G ot 1 -E) 



0. 



(18) 



To perform a type III approximation on F[E], we take 
advantage that it is universal, i.e., that it depends only 
on the interaction part of the Hamiltonian and not on 
the one-body part. This follows from the universal char- 
acter of its Legendre transform $[G]. We thus evaluate 
F[E] exactly for a Hamiltonian H' that shares the same 



interaction part as the Hubbard Hamiltonian, but that is 
exactly solvable. This Hamiltonian W is taken as a clus- 
ter decomposition of the original problem, i.e., we tile the 
infinite lattice into identical, disconnected clusters that 
can be solved exactly. Examples of such tilings are given 
in Fig. 13. Denoting the corresponding quantities with a 
prime, we obtain, 



O t ,[E'] = F[E'] - Trln(-GoJ + E'). 
from which we can extract -F[E']. It follows that 



(19) 



fi t [£'] = n t , [£'] + Tr ln(-G^ + E') - Tr ln(-G, 



ot 



-E'). 

The type III approximation comes from the fact that 
the self-energy £' is restricted to the exact self-energy of 
the cluster problem H' , so that variational parameters 
appear in the definition of the one-body part of H' . 

In practice, we look for values of the cluster one-body 
parameters t' such that <5Qt[E']/W = 0. It is useful 
for what follows to write the latter equation formally, 
although we do not use it in actual calculations. Given 
that fit'[E'] is the actual grand potential evaluated for 
the cluster, dflt'[S']/dt' is canceled by the explicit t' 
dependence of Trln(— G^J + £') and we are left with 

5fit[E'] <5E' 







<5£' 
-Tr 



St' 



1 



Got _ s ' 



8E_ 
St' 



(21) 



Given that the clusters corresponding to t' are discon- 
nected and that translation symmetry holds on the su- 
perlattice of clusters, each of which contains N c sites, the 
last equation may be written 



EE 



N 



1 



N c VG ot i-E'(z^W 




(St 1 



0. (22) 



3. Variational cluster perturbation theory, or variational 
cluster approximation 

In Variational Cluster Perturbation Theory (VCPT), 
more aptly named the Variational Cluster Approach 
(VCA), solutions to the Euler equations (22) are found by 
looking for numerical minima (or more generally, saddle- 
points) of the functional. Typically, the VCA cluster 
Hamiltonian H' will have the same form as H except 
that there is no hopping between clusters and that long- 
range order is allowed by adding some Weiss fields, for 
instance like in Eq. (37) below. The hopping terms and 
chemical potential within H' may also be treated like ad- 
ditional variational parameters. In contrast with Mean- 
Field theory, these Weiss fields are not mean fields, in 
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the sense that they do not coincide with the correspond- 
ing order parameters. The interaction part of H (or 
H') is not factorized in any way and short-range cor- 
relations arc treated exactly. In fact, the Hamiltonian 
H is not altered in any way; the Weiss fields are intro- 
duced to let the variational principle act on a space of 
self-energies that includes the possibility of specific long- 
range orders, without imposing those orders. Indeed, the 
more naturally an order arises in the system, the smaller 
the Weiss field needs to be, and one observes that the 
strength of the Weiss field at the stationary point of the 
self-energy functional generally decreases with increasing 
cluster size, as it should since in the thermodynamic limit 
no Weiss field should be necessary to establish order. 



4- Cellular dynamical mean-field theory 



number of degrees of freedom in the bath is insufficient. 
Instead, one adopts the following self-consistent scheme: 
(1) one starts with a guess value of the bath parameters 
(V^ a ,e a ) and solves the cluster Hamiltonian H' numeri- 
cally. (2) One then calculates the combination 



E 



+ E'(iw„) (26) 



and (3) minimizes the following canonically invariant dis- 
tance function: 



d= ^ (iu n + //-£'- t(iw n ) - G 1 ) 



(27) 



The Cellular dynamical mean-field theory (CDMFT) 
is obtained by including in the cluster Hamiltonian H' a 
bath of uncorrelated electrons that somehow must mimic 
the effect on the cluster of the rest of the lattice. Explic- 
itly, H' takes the form 

/i.,ZA(X fl 

(23) 

where a aa annihilates an electron of spin a on a bath 
orbital labelled a. The bath is characterized by the en- 
ergy of each orbital (e a ) and the bath-cluster hybridiza- 
tion matrix V^ a . This representation of the environment 
through an Anderson impurity model was introduced in 
Ref. 65 in the context of DMFT (i.e., a single site). The 
effect of the bath on the electron Green function is en- 
capsulated in the so-called hybridization function 

IV(u,) = £^^ (24) 

which enters the Green function as 

[G'-'U = " + M - C - - ( 25 ) 

Moreover, the CDMFT does not look for a strict solu- 
tion of the Euler equation (22), but tries instead to set 
each of the terms between brackets to zero separately. 
Since the Euler equation (22) can be seen as a scalar 
product, CDMFT requires that the modulus of one of the 
vectors vanish to make the scalar product vanish. From 
a heuristic point of view, it is as if each component of 
the Green function in the cluster were equal to the cor- 
responding component deduced from the lattice Green 
function. This clearly reduces to single site DMFT when 
there is only one lattice site. 

When the bath is discretized, i.e., is made of a finite 
number of bath "orbitals" , the left-hand side of Eq. (22) 
cannot vanish separately for each frequency, since the 



over the set of bath parameters (changing the bath pa- 
rameters at this step does not require a new solution of 
the Hamiltonian H', but merely a recalculation of the hy- 
bridization function T). The bath parameters obtained 
from this minimization are then put back into step (1) 
and the procedure is iterated until convergence. 

In practice, the distance function (27) can take 
various forms, for instance by adding a frequency- 
dependent weight in order to emphasize low-frequency 
properties 17 ' 59,66 or by using a sharp frequency cutoff. 67 
These weighting factors can be considered as rough ap- 
proximations for the missing factor ST,' (ioj n )/St' in the 
Euler equation (22). The frequencies are summed over on 
a discrete, regular grid along the imaginary axis, defined 
by some fictitious inverse temperature (3, typically of the 
order of 20 or 40 (in units of t^ 1 ). Even when the total 
number of cluster plus bath sites in CDMFT equals the 
number of sites in a VCA calculation, CDMFT is much 
faster than the VCA since the minimization of a grand 
potential functional requires many exact diagonalizations 
of the cluster Hamiltonian H' . 

The final lattice Green function from which one com- 
putes observable quantities may be obtained by periodiz- 
ing the self-energy, as in Ref. 58 or in the CPT manner 
described above in Eq. (11). We prefer the last approach 
because it corresponds to the Green function needed to 
obtain the density from dQ/d/j, = — Tr(G) and also be- 
cause periodization of the self-energy gives additional un- 
physical states in the Mott gap 68 (see also Ref. 60). 



5. The Dynamical cluster approximation 

The DCA 55 cannot be formulated within the self- 
energy functional approach. 159 It is based on the idea 
of discretizing irreducible quantities, such as the energy, 
in reciprocal space. It is believed to converge faster for 
q = quantities whereas CDMFT converges exponen- 
tially fast for local quantities. 69-71 
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FIG. 14: The spectral function of the U — » co limit of the 
one- dimensional Hubbard model, as calculated from (a) an 
exact diagonalization of the Hubbard model with U/t — 100 
on a periodic 12- site cluster; (b) the same, but with CPT, on 
a 12-site cluster with open boundary conditions; (c) the exact 
solution, taken from Ref. 72; beware: the axes are oriented 
differently. In (a) and (b) a finite width r] has been given to 
peaks that would otherwise be Dirac ^-functions. 



6. Benchmarks for quantum cluster approaches 

Since DMFT becomes exact in infinite dimension, the 
most difficult challenge for cluster extensions of this ap- 
proach is in one dimension. In addition, exact results to 
compare with exist only in one dimension so it is mostly 
in d = 1 that cluster methods have been checked. In 
d = 2 there have also been a few comparisons with QMC 
as we shall discuss. 

CPT has been checked 68 for example by comparing 
with exact results 72 for the spectral function at U — > oo 
in d = 1 as shown in Fig. 14. Fig. 15 shows the chemical 
potential as a function of density for various values of 
U . Fig. 16 shows the convergence rates for the total 
energy and for the double occupancy in the d = 1 half- 
filled model. Clearly, there is a dramatic improvement 
compared with exact diagonalizations. 




FIG. 15: Chemical potential as a function of density in the 
one-dimensional Hubbard model, as calculated by CPT (from 
Ref. 62). The exact, Bethe-Ansatz result is shown as a solid 
line. 
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FIG. 16: Top: Comparison (expressed in relative difference) 
between the ground-state energy density of the half-filled, 
one- dimensional Hubbard model calculated from the exact, 
Bethe-Ansatz result. The results are displayed as a function 
of the hopping t, for U = 2t and various cluster sizes L (con- 
nected symbols). For comparison, the exact diagonalization 
values of finite clusters with periodic boundary conditions are 
also shown (dashed lines) for L — 8 and L = 12. Bottom: 
Same for the double occupancy. An extrapolation of the re- 
sults to infinite cluster size (L — > oo) using a quadratic fit in 
terms of 1/L is also shown, and is accurate to within 0.5%. 
Taken from Ref. 62. 



The main weakness of CPT is that it cannot take 
into account tendency towards long-range order. This 
is remedied by VCPT, as shown in Fig. 17 where CPT, 
VCPT are both compared with QMC as a benchmark. 
Despite this agreement, we should stress that long wave 
length fluctuations are clearly absent from cluster ap- 
proaches. Hence, the antiferromagnctic order parameter 
at half-filling, for example, does not contain the effect of 
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FIG. 17: Ground state energy of the half-filled, two- 
dimensional Hubbard model (t = 1) as a function of U, as 
obtained from various methods: Exact diagonalization (ED), 
CPT and VCPT on a 10-site cluster, quantum Monte Carlo 
(QMC) and variational Monte Carlo (VMC). Taken from 
Ref. 73. 




FIG. 18: CDMFT calculation on a 2 x 2 cluster with 8 bath 
sites of the density as a function of the chemical potential 
in the one- dimensional Hubbard model for U = 4i, as com- 
pared with the exact solution, DMFT and other approxima- 
tion schemes. Taken from Ref. 74. 



temperature defining the Matsubara frequencies to sum 
over. The final results are not completely insensitive to 
the choice of fictitious temperature or frequency- weighing 
scheme but are usually considered reliable and consistent 
with each other when (3 lies between 20 and 40. The cut- 
off procedures have been discussed in Ref. 67. 

The relative merits of DCA and CDMFT have been 
discussed for example in Refs 69-71,75,76. Briefly speak- 
ing, convergence seems faster in DCA for long wave 
length quantities but CDMFT is faster (exponentially) 
for local quantities. 



III. RESULTS AND CONCORDANCE 
BETWEEN DIFFERENT METHODS 

In this section, we outline the main results we obtained 
concerning the pseudogap and d-wave superconductivity 
in the two-dimensional Hubbard model. Quantum clus- 
ter approaches are better at strong coupling while TPSC 
is best at weak coupling. Nevertheless, all these methods 
are non-perturbative, the intermediate coupling regime 
presenting the physically most interesting case. But it 
is also the regime where we have the least control over 
the approximations. As we will show, it is quite satis- 
fying that, at intermediate coupling, weak-coupling and 
strong-coupling approaches give results that are nearly 
in quantitative agreement with each other. This gives 
us great confidence into the validity of the results. As 
an example of concordance, consider the fact that to ob- 
tain spectral weight near (jr/2, n/2) at optimal doping in 
the electron-doped systems, U has to be roughly 6t. For 
larger U, (U — 8t in CPT) that weight, present in the 
experiments, disappears. Smaller values of U (U = 4t 
in CPT) do not lead to a pseudogap. Other examples 
of concordance include the value of the superconducting 
transition temperature T c obtained with DCA and with 
TPSC as well as the temperature dependence of souble 
occupancy obtained with the same two methods. 



zero-point long wave length transverse spin fluctuations. 
This is discussed for example in the context of Fig. 9 of 
Ref. 73. 

CDMFT corrects the difficulties of CPT near half- 
filling by reproducing the infinite compressibility pre- 
dicted by the Bethe ansatz in one dimension as shown 
in Fig. 18. 74 Detailed comparisons between the local and 
near-neighbor Green functions 66 ' 74 have been performed. 
One should note that these results, obtained from ex- 
act diagonalization, also need the definition of a distance 
function (See Eq. (27) above) that helps find the best 
bath parametrization to satisfy the self-consistency con- 
dition. This measure forces one to define calculational 
parameters such as a frequency cutoff and an fictitious 



A. Weak and strong-coupling pseudogap 

To understand the pseudogap it is most interesting to 
consider both hole and electron-doped cuprates at once. 
This means that we have to include particle-hole sym- 
metry breaking hoppings, t' and t" . We will see in the 
present section that it is possible to obtain a pseudogap at 
strong coupling without a large correlation length in the 
particle- hole or in the particle-particle channels. By con- 
trast, at weak coupling one does need a long-correlation 
length and low dimension. So there appears to be theo- 
retically two different mechanisms for the pseudogap. 
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FIG. 19: Single particle spectral weight, as a function of en- 
ergy u> in units of t, for wave vectors along the high-symmetry 
directions shown in the inset, (a) CPT calculations on a 3 x 4 
cluster with ten electrons (17% hole doped), (b) the same 
as (a), with 14 electrons (17% electron doped). In all cases 
t' = — 0.3t and t" = 0.2t. A Lorentzian broadening 77 = 0.12t 
is used to reveal the otherwise delta peaks. From Ref. 77. 




FIG. 20: Onset of the pseudogap as a function of U corre- 
sponding to Fig. 19, taken from Ref. 77. Hole-doped case on 
the left, electron-doped case on the right panel 



1, CPT 

The top panel in Fig. 19 presents the single-particle 
spectral weight, A(k, ui) or imaginary part of the single- 
particle Green function, for the model with t' — — 0.3t, 
t" = Q.2t in the hole-doped case, for about 17% doping. 77 




FIG. 21: MDC from CPT in the t-t' = -Q.3t, t" = 0.2t 
Hubbard model, taken from Ref. 77. 



Each curve is for a different wave vector (on a trajectory 
shown in the inset) and is plotted as a function of fre- 
quency in units of t. This kind of plot is known as En- 
ergy Dispersion Curves (EDC). It is important to point 
out that the theoretical results are obtained by broaden- 
ing a set of delta function, so that the energy resolution 
is r\ = 0.12t corresponding roughly to the experimental 
resolution we will compare with in the next section. At 
small U = 2t on the top panel of Fig. 19, one recovers a 
Fermi liquid. At large U, say U = 8t, the Mott gap at 
positive energy is a prominent feature. The pseudogap 
is a different feature located around the Fermi energy. 
To see it better, we present on the left-hand panel of 
Fig. 20 a blow-up in the vicinity of the Fermi surface 
crossing occurring near (tt, 0). Clearly, there is a mini- 
mum in AQt, ui) at the Fermi-surface crossing when U is 
large enough instead of a maximum like in Fermi liquid 
theory. 

It is also possible to plot A(k, ui) at fixed frequency for 
various momenta. They are so-called Momentum Dis- 
persion Curves (MDC). In Fig. 21 we take the Fermi 
energy uo = 0, and we plot the magnitude of the single- 
particle spectral weight in the first quadrant of the Bril- 
louin zone using red for high-intensity and blue for low 
intensity. The figure shows that, in the hole-doped case 
(top panel), weight near (tt/2, tt/2) survives while it tends 
to disappear near (tt, 0) and (0, tt). That pseudogap phe- 
nomenon is due not only to large U but also to the fact 
that the line that can be drawn between the points (tt, 0) 
and (0, tt) crosses the Fermi surface. When there is no 
such crossing, one recovers a Fermi surface (not shown 
here). The (71-, 0) to (0, tt) line has a double role. It is the 
antiferromagnetic zone boundary, as well as the line that 
indicates where umklapp processes become possible, i.e., 
the line where we can scatter a pair of particles on one 
side of the Fermi surface to the other side with loss or 
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FIG. 22: Right: The corresponding EDC in the t-t'-t" Hub- 
bard model, calculated on a 16-site cluster in CPT, at n = 
9/8. Inset: the pseudogap. Left: The corresponding Fermi 
energy momentum distribution curve. 
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FIG. 23: EDC in the t-t'-t" Hubbard model, with t' = -0.3* 
and t ' = 0.2t, calculated on a 8-site cluster for U = 8t 
in VCPT. d-wave superconductivity is present in the holde- 
doped case (left) and both antiferromagnetism and d-wave 
superconductivity in the electron-doped case. The resolution 
is not large enough in the latter case to see the superconduct- 
ing gap. The Lorentzian broadening is 0.2t. From Ref. 19. 
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FIG. 24: MDC in the t-t' , U = 8t Hubbard model, calculated 
on a 4-site cluster in CDMFT. Energy resolution, r\ — O.lt 
(left and middle). Left: Hole-doped dSC (t' = -0.3t,n = 
0.96), Middle: Electron-doped dSC (t' = 0.3t,n = 0.93), 
Right: Same as middle with r\ = 0.02t. Note the particle-hole 
transformation in the electron-doped case. From Ref. 17. 




FIG. 25: EDC in the t-t 1 = 0, U = 8t Hubbard model, cal- 
culated on a 4-site cluster in CDMFT. Top: normal (para- 
magnetic) state for various densities. Bottom: same for the 
antiferromagnetic state. From Ref. 67. 



gain of a reciprocal lattice vector. Large scattering rates 
explain the disappearance of the Fermi surface. 77 We also 
note that the size of the pseudogap in CPT, defined as 
the distance between the two peaks, does not scale like 
J = 4t 2 /U at large coupling. It seems to be very weakly 
U dependent, its size being related to t instead. This 
result is corroborated by CDMFT. 67 

The EDC for the electron-doped case is shown on the 
bottom panels of Fig. 19 near optimal doping again. This 
time, the Mott gap appears below the Fermi surface so 
that the lower Hubbard band becomes accessible to ex- 
periment. The EDC in Fig. 22 shows very well both the 
Mott gap and the pseudogap. Details of that pseudogap 
can be seen both in the inset of Fig. 22 or on the right- 
hand panel of Fig. 20. While in the hole-doped case the 
MDC appeared to evolve continuously as we increase U 
(top panel of Fig. 21), in the electron-doped case (bot- 
tom panel) the weight initially present near (7r/2, 7t/2) at 
U = At disappears by the time we reach U — 8t. 

In Fig. 23 we show, with the same resolution as CPT, 
the MDC for VCPT. 19 In this case the effect of long-range 
order is included and visible but, at this resolution, the 
results are not too different from those obtained from 
CPT in Fig. 21. 



2. CDMFT and DCA 

CDMFT 17 gives MDC that, at comparable resolution, 
i] = O.lt, are again compatible with CPT and with 
VCPT. The middle panel in Fig. 24 is for the electron- 
doped case but with a particle-hole transformation so 
that t' = +0.3t and k — > k + (tt,tt). Since there is a 
non-zero d-wave order parameter in this calculation, im- 
proving the resolution to r\ = 0.02£ reveals the d-wave 
gap, as seen in the right most figure. 

It has been argued for a while in DCA that there 
is a mechanism whereby short-range correlations at 
strong coupling can be the source of the pseudogap 
phenomenon. 78 To illustrate this mechanism in CDMFT, 
we take the case t' = t" = and compare in Fig. 25 
the EDC for U = 8t without long-range order (top pan- 
els) and with long-range antiferromagnetic order (bot- 
tom panels). 67 The four bands appearing in Figs 25a and 
25d are in agreement with what has been shown 73,79 ' 80 
with CPT, VCPT and QMC in Fig. 26. Evidently there 
are additional symmetries in the antiferromagnetic case. 
The bands that are most affected by the long-range or- 
der are those that are closest to the Fermi energy, hence 
they reflect spin correlations, while the bands far from 
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FIG. 26: EDO in the Hubbard model, U = 8t, t' 
lated in CPT, VCPT and QMC. From Ref. 73. 
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the Fermi energy seem less sensitive to the presence of 
long-range order. These far away bands are what is left 
from the atomic limit where we have two dispersionless 
bands. As we dope, the chemical potential moves into 
the lower band closest to the Fermi energy. When there 
is no long-range order (Figs 25b and 25c) the lower band 
closest to the Fermi energy moves very close to it, at the 
same time as the upper band closest to the Fermi energy 
looses weight, part of it moving closer to the Fermi en- 
ergy. These two bands leave a pseudogap at the Fermi 
energy 82,83 , although we cannot exclude that increasing 
the resolution would reveal a Fermi liquid at a very small 
energy scale. In the case when there is long-range anti- 
ferromagnetic order, (Figs 25e and 25f) the upper band 
is less affected while the chemical potential moves in the 
lower band closest to the Fermi energy but without creat- 
ing a pseudogap, as if we were doping an itinerant antifer- 
romagnet. It seems that forcing the spin correlations to 
remain short range leads to the pseudogap phenomenon 
in this case. When a pseudogap appears, it is created 
again by very large scattering rates. 67 



3. TPSC (including analytical results) 

In Hartree-Fock theory, double occupancy is given by 
n 2 /4 and is independent of temperature. The correct 
result does depend on temperature. One can observe 
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FIG. 27: Double occupancy (n^ni), in the two-dimensional 
Hubbard model for n = 1, as calculated from TPSC (lines 
with x) and from DCA (symbols) from Ref. 81. Taken from 
Ref. 84. 
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FIG. 28: MDC at the Fermi energy for the two-dimensional 
Hubbard model for U = 6.25, t' = -0.175i,t" = 0.05t at 
various hole dopings, as obtained from TPSC. The far left 
from Ref. 85 is the Fermi surface plot for 10% hole-doped 
Caa-^Na.CuOaCla. 



in Fig. 27 the concordance between the results for the 
temperature-dependent double occupancy obtained with 
DCA and with TPSC 84 for the t' = t" = Q model. We 
have also done extensive comparisons between straight 
QMC calculations and TPSC. 38 The downturn at low 
temperature has been confirmed by the QMC calcula- 
tions. It comes from the opening of the pseudogap due 
to antiferromagnetic fluctuations, as we will describe be- 
low. The concomitant increase in the local moment corre- 
sponds to the decrease in double-occupancy. There seems 
to be a disagreement at low temperature between TPSC 
and DCA at U = 2t. In fact TPSC is closer to the di- 
rect QMC calculation. Since we expect quantum cluster 
methods in general and DCA in particular to be less ac- 
curate at weak coupling, this is not too worrisome. At 
U = At the density of states obtained with TPSC and 
with DCA at various temperatures are very close to each 
other. 84 We stress that as we go to temperatures well 
below the pseudogap, TPSC becomes less and less ac- 
curate, generally overemphasizing the downfall in double 
occupancy. 

We will come back to more details on the predictions 
of TPSC for the pseudogap, but to illustrate the concor- 
dance with quantum cluster results shown in the previ- 
ous subsection, we show in Fig. 28 MDC obtained at the 
Fermi energy in the hole doped case for t' = — 0.175i, 
t" = 0.05t. Again there is quasi-particle weight near 
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FIG. 29: MDC at the Fermi energy in the electron-doped case 
with t' = -0.1754, t" = 0.05t and two different U's, U = 6.254 
and U = 5.754 obtained from TPSC. The first column is the 
corresponding experimental plots at 10% and 15% doping in 
Ref. 86. From Refs. 40,87. 
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FIG. 30: Cartoon explanation of the pseudogap in the weak- 
coupling limit. Below the dashed crossover line to the renor- 
malized classical regime, when the antiferromagnetic corre- 
lation length becomes larger than the thermal de Broglie 
wave length, there appears precursors of the zero-temperature 
Boboliubov quasiparticles for the long-range ordered antifer- 
romagnet. 



(7r/2,7r/2) and a pseudogap near (it, 0) and (0, ir). How- 
ever, as we will discuss below, the antiferromagnetic cor- 
relation length necessary to obtain that pseudogap is too 
large compared with experiment. The electron-doped 
case is shown in Fig. 29 near optimal doping and for 
different values of U. As U increases, the weight near 
(■7r/2,7r/2) disappears. That is in concordance with the 
results of CPT shown in Fig. 21 where the weight at 
that location exists only for small U. That also agrees 
with slave-boson calculations 88 that found such weight 
for U — 6t and it agrees also with one-loop calculations 89 
starting from a Hartree-Fock antiferromagnetic state that 
did not find weight at that location for U — 8t. The 
simplest Hartree-Fock approach 21 ' 90 yields weight near 
{it/2, 7t/2) only for unreasonably small values of U . 
A cartoon explanation of the pseudogap is given in 



Fig. 30. At high temperature we have a Fermi liquid, 
as illustrated in panel I. Now, suppose we start from a 
ground state with long-range order as in panel III, in 
other words at a filling between half- filling and n c . In 
the Hartree-Fock approximation we have a gap and the 
fcrmion creation-annihilation operators now project on 
Bogoliubov- Valentin quasiparticles that have weight at 
both positive and negative energies. In two dimensions, 
the Mermin- Wagner theorem means that as soon as we 
raise the temperature above zero, long-range order dis- 
appears, but the antiferromagnetic correlation length £ 
remains large so we obtain the situation illustrated in 
panel II, as long as £ is much larger than the thermal de 
Broglie wave length fy, = vf/(ttT) in our usual units. At 
the crossover temperature Tx then the relative size of £ 
and ^th changes and we recover the Fermi liquid. We now 
proceed to sketch analytically where these results come 
from starting from finite temperature. Details and more 
complete formulae may be found in Refs. 23,24,33,37. 
Note also that a study starting from zero temperature 
has also been performed in Ref. 91. 

First we show how TPSC recovers the Mermin- Wagner 
theorem. Consider the self-consistency conditions given 
by the local moment sum rule Eq. (4) together with the 
expression for the spin-susceptibility Eq. (2) and U sp in 
Eq. (5). First, it is clear that if the left-hand side of the 
local moment sum rule Eq. (4) wants to increase because 
of proximity to a phase transition, the right-hand side can 
do so only by decreasing (n-\n\) which in turns decreases 
U sp through Eq. (5) and moves the system away from the 
phase transition. This argument needs to be made more 
precise to include the effect of dimension. First, using 
the spectral representation one can show that every term 
of Xsp(q, iQn) is positive. Near a phase transition, the 
zero Matsubara frequency component of the susceptibil- 
ity begins to diverge. On can check from the real-time 
formalism that the zero-Matsubara frequency contribu- 
tion dominates when the characteristic spin fluctuation 
frequency uj sp ~ £ -2 becomes less than temperature, the 
so-called renormalized-classical regime. We isolate this 
contribution on the left-hand side of the local moment 
sum rule and we move the contributions from the non- 
zero Matsubara frequencies, that are non-divergent, on 
the right-hand side. Then, converting the wave vector 
sum to an integral and expanding the denominator of the 
susceptibility around the wave vector where the instabil- 
ity would occurs to obtain an Ornstein-Zernicke form, 
the local moment sum rule Eq. (4) can be written in the 
form 



T 



l dq- 



1 



r 



= C(T). 



(28) 



The constant on the right-hand side contains only non- 
singular contributions and £ -2 contains U sp that we want 
to find. From the above equation, one finds immediately 
that in d = 2, £ w exp(C(T)/T) so that the correla- 
tion length diverges only at T — 0. In three dimensions, 
isotropic or not, exponents correspond to those of the 
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FIG. 31: MDC plots at the Fermi energy (upper) and corre- 
sponding scattering rates (lower) obtained from TPSC. The 
red lines on the upper panel indicate the region where the 
scattering rate in the corresponding lower panels is large. 



N = oo universality class. 46 

To see how the pseudogap opens up in the single- 
particle spectral weight, consider the expression (6) for 
the self-energy. Normally one has to do the sum over 
bosonic Matsubara frequencies first, but the zero Mat- 
subara frequency contribution has the correct asymp- 
totic behavior in fermionic frequencies ik n so that one 
can once more isolate on the right-hand side the zero 
Matsubara frequency contribution. This is confirmed by 
the real-time formalism 23 (See also Eq. (36) below). In 
the renormalized classical regime then, we have 



E(kjr, ifcn) oc T J q' 1 - 1 



dq- 



1 



1 



£ 2 ik n - £k F +Q+ q 



(29) 



where Q is the wave vector of the instability. Hence, 
when £k F +Q = 0, in other words at hot spots, we find af- 
ter analytical continuation and dimensional analysis that 



ImE fl (k F ,0) oc -ttT J < 



jd-l 



q±dqn- T - 
q± 



l 



ttT 



3-d 



Z2 S ( V FQ\\) 

(30) 
(31) 



Clearly, in d = 4, ImE fi (kF,0) vanishes as temperature 
decreases, d = 3 is the marginal dimension and in d = 2 
we have that ImE^kp ,0) oc £/£th that diverges at zero 
temperature. In a Fermi liquid that quantity vanishes at 
zero temperature. A diverging Im S fi (kj?, 0) corresponds 
to a vanishingly small A(hp, uj = 0) as we can see from 



A(k,uj) 



-2ImE fl (k F ,w) 



{uj - e k - Re E fl (k F) uj)) 2 + Im E fl (k F , uj) 2 ' 

(32) 

To see graphically this relationship between the loca- 
tion of the pseudogap and large scattering rates at the 



Fermi surface, we draw in Fig. 31 both the Fermi surface 
MDC and, in the lower panels, the corresponding plots 
for ImE fl (k, 0). Note that at stronger [7 the scattering 
rate is large over a broader region, leading to a depletion 
of A(k,u>) over a broader range of k values. 

An argument for the splitting in two peaks seen in 
Figs. 6 and 30 is as follows. Consider the singular renor- 
malized contribution coming from the spin fluctuations in 
Eq. (29) at frequencies w ^> Wf£ _1 - Taking into account 
that contributions to the integral come mostly from a 
region q < this expression leads to 



ReE fl (k F ,u;) = (t J q^dq^-^ 



1 



A 2 



w - £k F +Q 



£ 2 J ik n - £k F +Q 

(33) 



which, when substituted in the expression for the spectral 
weight (32) leads to large contributions when 



w - £k - 



A 



W _ £ k F +Q 



= 



(34) 



or, equivalently, 



(£k + £k F+Q ) ± Vgk - £k F+Q ) 2 + 4A 2 
w = 2 ' ' ) 

which corresponds to the position of the hot spots in 
Fig. 29 for example. 

Note that analogous arguments hold for any fluctu- 
ation that becomes soft, 23 including superconducting 
ones. 41,43 The wave vector Q would be different in each 
case. 



4- Weak- and strong- coupling pseu 



The CPT results of Figs. 19 and 22 clearly show that 
the pseudogap is different from the Mott gap. At finite 
doping, the Mott gap remains a local phenomenon, in 
the sense that there is a region in frequency space that is 
not tied to uj = where for all wave vectors there are no 
states. The peudogap by contrast is tied to uj = and 
occurs in regions nearly connected by (tt, 7t), whether we 
are talking about hole- or about electron-doped cuprates. 
That the phenomenon is caused by short-range correla- 
tions can be seen in CPT from the fact that the pseu- 
dogap is independent of cluster shape and size (most of 
the results were presented for 3x4 clusters and we did 
not go below size 2x2). The antiferromagnetic cor- 
relations and any other two-particle correlations do not 
extend beyond the size of the lattice in CPT. Hence, the 
pseudogap phenomenon cannot be caused by antiferro- 
magnetic long-range order since no such order exists in 
CPT. This is also vividly illustrated by the CDMFT re- 
sults in Fig. 25 that contrast the case with and without 
antiferromagnetic long-range order. The CDMFT results 
also suggest that the pseudogap appears in the bands 
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that are most affected by antifcrromagnctic correlations 
hence it seems natural to associate it with short-range 
spin correlations. The value of t' has an effect, but it 
mostly through the fact that it has a strong influence on 
the relative location of the antiferromagnetic zone bound- 
ary and the Fermi surface, a crucial factor determining 
where the pseudogap is. All of this as well as many re- 
sults obtained earlier by DCA 78 suggest that there is a 
strong coupling mechanism that leads to a pseudogap in 
the presence of only short-range two-body correlations. 
However, the range cannot be zero. Only the Mott gap 
appears at zero range, thus the pseudogap is absent in 
single-site DMFT. 

In the presence of a pseudogap at strong coupling 
(U > 8t), wave vector is not, so to speak, such a bad 
quantum number in certain directions. In other words 
the wave description is better in those directions. In 
other directions that are "pseudogapped" , it is as if the 
localized, or particle description was better. This com- 
petition between wave and particle behavior is inherent 
to the Hubbard model. At the Fermi surface in low di- 
mension, it seems that this competition is resolved by 
dividing (it is a crossover, not a real division) the Fermi 
surface in different sections. 

There is also a weak-coupling mechanism for the pseu- 
dogap. This has been discussed at length just in the pre- 
vious section on TPSC. Another way to rephrase the cal- 
culations of the previous section is in the real frequency 
formalism. There one finds 23 that 

£"^(k F ,c) 

K / W^ 1 1 v [n(w ' ) + /(w + w ' )]x - (q;w,) 

(36) 

so that if the characteristic spin fluctuation frequency 
in x" p (q;a/) is much larger than temperature, then 
[n{oj') + f(u)+u>')] can be considered to act like a window 
of size lo or T and x' S p(q; w ') can be replaced by a func- 
tion of q times lo' which immediately leads to the Fermi 
liquid result [lo 2 + (irT) 2 ]. In the opposite limit where the 
characteristic spin fluctuation frequency in x" p (q;o/) is 
much less than temperature, then it acts as a window 
narrower than temperature and [n{io') + f(to + lo')] can 
be approximated by the low frequency limit of the Bose 
factor, namely T/lo'. Using the thermodynamic sum rule, 
that immediately leads to the result discussed before in 
Eq.(31), ImS(k F ,0) oc (TrT/v' F )f- d . This mechanism 
for the pseudogap needs long correlation lengths. In 
CPT, this manifests itself by the fact that the apparent 
pseudogap in Fig. 21 at U — At is in fact mostly a depres- 
sion in spectral weight that depends on system size and 
shape. In addition, in contrast to the short-range strong- 
coupling mechanism, at weak coupling the pseudogap is 
more closely associated with the intersection of the anti- 
ferromagnetic zone boundary with the Fermi surface. 

Which mechanism is important for the cuprates will 
be discussed below in the section on comparisons with 
experiments. 



B. d-wave superconductivity 

The existence of d-wave superconductivity at weak 
coupling in the Hubbard model mediated by the 
exchange of antiferromagnetic fluctuations 92 ' 93 had 
been proposed even before the discovery of high- 
temperature superconductivity. 6-8 At strong-coupling, 
early papers 13,1 also proposed that the superconductiv- 
ity would be d-wave. The issue became extremely contro- 
versial, and even recently papers have been published 12 
that suggest that there is no d-wave superconductivity 
in the Hubbard model. That problem could have been 
solved very long ago through QMC calculations if it had 
been possible to do them at low enough temperature. 
Unfortunately, the sign problem hindered these simula- 
tions, and the high temperature results 5,44 ' 95,96 were not 
encouraging: the d-wave susceptibility was smaller than 
for the non-interacting case. Since that time, numerical 
results from variational QMC, 4,97 exact diagonalization 16 
and other numerical approaches 98 for example, suggest 
that there is indeed d-wave superconductivity in the Hub- 
bard model. 

In the first subsection, we show that VCPT leads to 
a zero-temperature phase diagram for both hole and 
electron-doped systems that does show the basic features 
of the cuprate phase diagram, namely an antiferromag- 
netic phase and a d-wave superconducting phase in dop- 
ing ranges that are quite close to experiment 19 (The fol- 
lowing section will treat in more detail comparisons with 
experiment). The results are consistent with CDMFT. 17 
The fall in the d-wave superconducting order parameter 
near half-filling is associated with the Mott phenomenon. 
The next subsection stresses the instability towards d- 
wave superconductivity as seen from the normal state 
and mostly at weak coupling. We show that TPSC can 
reproduce available QMC results and that its extrapo- 
lation to lower temperature shows d-wave superconduc- 
tivity in the Hubbard model. The transition tempera- 
ture found at optimal doping 40 for U — At agrees with 
that found by DCA, 18 a result that could be fortuitous. 
But again the concordance between theoretical results 
obtained at intermediate coupling with methods that are 
best at opposite ends of the range of coupling strengths 
is encouraging. 

1. Zero-temperature phase diagram 

In VCPT, one adds to the cluster Hamiltonian the 
terms 19 

H' M = MY j e iQ ^(n^-n tli ) (37) 
H' D =Dj29„Ac^c vl +R.c.) (38) 

with M and D are respectively antiferromagnetic and d- 
wave superconducting Weiss fields that are determined 
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FIG. 32: Antiferromagnetic (bottom) and d-wave (top) order 
parameters for U = 8t, t' — — 0.3t t" — 0.2t as a function of 
the electron density (n) for 2 x 3, 2 x 4 and 10-site clusters, 
calculated in VCPT. Vertical lines indicate the first doping 
where only d-wave order is non- vanishing. From Ref. 19. 
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FIG. 33: d-wave order parameter as a function of n for various 
values of t' , calculated in CDMFT on a 2 x 2 cluster for U = 
St. The positive t' case corresponds to the electron-doped 
case when a particle-hole transformation is performed. From 
Ref. 17. 



self-consistently and g^v equal to ±1 on near-neighbor 
sites following the d-wave pattern. We recall that the 
cluster Hamiltonian should be understood in a varia- 
tional sense. Fig. 32 summarizes, for various cluster sizes, 
the results for the d-wave order parameter D and for the 
antiferromagnetic order parameter Mq for a fixed value of 
U = St and the usual hopping parameters t' = — 0.3t and 
t" = 0.2t. The results for antiferromagnetism are quite 
robust and extend over ranges of dopings that correspond 
quite closely to those observed experimentally. Despite 
the fact that the results for d-wave superconductivity still 
show some size dependence, it is clear that supercon- 
ductivity alone without coexistence extends over a much 
broader range of dopings on the hole-doped than on the 
electron-doped side as observed experimentally. VCPT 
calculations on smaller system sizes" but that include, 
for thermodynamic consistency, the cluster chemical po- 
tential as a variational parameter show superconductivity 
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FIG. 34: Antiferromagnetic (bottom) and d-wave (top) order 
parameters as a function of the electron density (n) for t' = 
— 0.3t t" = 0.2t and various values of U on a 8-site cluster, 
calculated in VCPT. From Ref. 19. 
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FIG. 35: d-wave order parameter as a function of n for various 
values of U, and t' = t" = calculated in CDMFT on a 4-site 
cluster. From Ref. 17. 



that extends over a much broader range of dopings. Also, 
for small 2x2 clusters, VCPT has stronger order parame- 
ter on the electron than on the hole-doped side, contrary 
to the results for the largest system sizes in Fig. 32. This 
is also what is found in CDMFT as can be seen in Fig. 33. 
It is quite likely that the zero-temperature Cooper pair 
size is larger than two sites, so we consider the results for 
2x2 systems only for their qualitative value. 

Concerning the question of coexistence with antifer- 
romagnetism, one can see that it is quite robust on the 
electron-doped side whereas on the hole-doped side, it 
is size dependent. That suggests that one should also 
look at inhomogeneous solutions on the hole-doped side 
since stripes are generally found experimentally near the 
regions where antiferromagnetism and superconductivity 
meet. 

Fig. 34 shows clearly that at strong coupling the size 
of the order parameter seems to scale with J, in other 
words it decreases with increasing U. This is also found 
in CDMFT, 17 as shown in Fig. 35 for t' = t" = 0. 

If we keep the antiferromagnetic order parameter to 
zero, one can check with both VCPT and CDMFT 
(Fig. 33) that the d-wave superconducting order parame- 
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FIG. 36: VCPT calculations for U = At, t' = t" = near 
half-filling on 2 x 4 lattice. Contrary to the strong coupling 
case, the d-wave order parameter Do survives all the way to 
half-filling at weak coupling, unless we also allow for antifer- 
romagnetism. 



ter goes to zero at half-filling. This is clearly due to Mott 
localization. Indeed, at smaller U — At for example, the 
order parameter docs not vanish at half-filling if we do 
not allow for long-range antiferromagnetic order, as il- 
lustrated in Fig. 35 for CDMFT. 17 The same result was 
found in VCPT, as shown in Fig. 36. 19 Restoring long- 
range antiferromagnetic order does however make the d- 
wavc order parameter vanish at half-filling. 

There are thus two ways to make d-wave superconduc- 
tivity disappear at half-filling, cither through long anti- 
ferromagnetic correlation lengths 94 or through the Mott 
phenomenon. In the real systems, that are Mott insula- 
tors and also antiferromagnets at half-filling, both effects 
can contribute. 



2. Instability of the normal phase 

In the introduction to this section, we alluded to QMC 
calculations for the d-wave susceptibility 5 ' 44,95 ' 96 Recent 
results 40 for that quantity as a function of doping for var- 
ious temperatures and for U = At, t' — t" = are shown 
by symbols in Fig. 8. For lower temperatures, the sign 
problem prevents the calculation near half-filling. Yet, 
the lowest temperature is low enough that a dome shape 
begins to appear. Nevertheless, comparison with the non- 
interacting case, shown by the top continuous line, leads 
one to believe that interactions only suppress d-wave su- 
perconductivity. We can easily understand why this is 
so. As we already know, the TPSC results obtained from 
Eq. (7) are very close to the QMC calculations, as shown 
by the solid lines in Fig. 8. In the temperature range 
of interest, the main contribution comes from the first 
term in Eq. (7). That term represents the contribution 
to the susceptibility that comes from dressed quasiparti- 
cles that do not interact with each other. Since dressed 
quasiparticles have a lifetime, a pair breaking effect, it is 
normal that this contribution to the interacting suscep- 
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FIG. 37: T c as a function of doping, 8 = 1 — n, for t = t = 
calculated in TPSC using the Thouless criterion. From 
Ref. 40. 



tibility leads to a smaller contribution than in the non- 
interacting case. At the lowest temperature, j3 = i/t, 
the vertex contribution represented by the second term 
in Eq. (7) accounts for about 20% of the total. It goes 
in the direction of increasing the susceptibility. As we 
decrease the temperature further in TPSC, that term be- 
comes comparable with the first one. Since the vertex in 
Eq. (7) accounts for the exchange of a single spin wave, 
equality with the first term signals the divergence of the 
series, as in 1/(1 — x) ~ 1 + x. The divergence of that 
series represents physically the instability of the normal 
phase to a d-wave superconducting phase. This is analo- 
gous to the Thouless criterion and hence it gives an up- 
per bound to T c . In other words, Berezinskii-Kostcrlitz- 
Thouless physics is not included. 

Fig. 37 shows the TPSC transition temperature for 
U = At and for U — 6t. As we move towards half- 
filling, located to the left of the diagram, starting from 
large dopings, T c goes up because of the increase in anti- 
ferromagnetic fluctuations. Eventually, T c turns around 
and decreases because of the opening of the pseudogap. 
Physically, when the pseudogap opens up, weight is re- 
moved from the Fermi level, the density of states be- 
comes very small, and pairing cannot occur any more. 
In the FLEX approximation 22 ' 100 that does not exhibit a 
pseudogap, 101 that downturn, observed already in QMC 
at high temperature, does not occur. We have observed 
that in cases where t' ^ so that the pseudogap opens 
only in a limited region around hot spots, the downturn 
can become less pronounced. 

The case n = 0.9 = 1 — 6 that corresponds to opti- 
mal doping for U = 4t in Fig. 37 has been studied by 
DCA. In an extensive and systematic study of the size 
dependence, Maier et al. ls established the existence of 
a d-wave instability at a temperature that coincides to 
within a few percent with the result in Fig. 37. Since very 
few vortices can fit within even the largest cluster sizes 
studied in Ref. 18, it is quite likely that the T c that they 
find does not include Berezinskii-Kosterlitz-Thouless ef- 
fects, just like ours. Despite the fact that, again, the 



21 



concordance between weak and strong coupling methods 
at intermediate coupling comforts us, the uncertainties in 
the results found with TPSC and DCA force us to also 
allow for a fortuitous coincidence. 



IV. COMPARISONS WITH EXPERIMENT 

The reduction of the real problem of high-temperature 
superconducting materials to a one-band Hubbard is a 
non-trivial one. It has been discussed already in the 
early days of high T c superconductivity. The notion of 
a Zhang-Rice singlet 102 emerged for hole doped systems. 
The mapping to a one-band model has been discussed 
in many references, 103 ' 104 and we do not wish to discuss 
this point further here. In fact it is far from obvious 
that this mapping is possible. It is known that about 0.5 
eV below the Fermi surface, that mapping fails in hole- 
doped systems. 10 Nevertheless, the one-band Hubbard 
model is in itself a hard enough problem for us. So it is 
satisfying to see that, in the end, it gives a picture that 
agrees with experiment in a quite detailed manner for the 
ARPES spectrum near the Fermi surface, for the phase 
diagram as well as for neutron scattering in cases where 
it can be calculated. 

Although wc will not come back on this point at all, we 
briefly mention that fitting the spin wave spectrum 105 for 
all energies and wave vectors at half filling in La2Cu04 
gives values of U,t,t',t" that are close to those used in 
the rest of this paper. 106-109 It is in this context that ring 
exchange terms are usually discussed. 

A. ARPES spectrum, an overview 

ARPES experiments have played a central role in the 
field of high-temperature superconductivity. We cannot 
expect to be able to present the vast experimental lit- 
erature on the subject. We refer the reader to a very 
exhaustive review 110 and to some less complete but re- 
cent ones. 111 ' 112 The main facts about ARPES have been 
summarized in Ref. 110. We comment on their main 
points one by one, using italics for our paraphrase of the 
reported experimental observations. 

(i) The importance of Mott Physics and the renormal- 
ization of the bandwidth from t to J for the undoped 
parent compounds. This renormalization was clear al- 
ready in early QMC calculations. 80,113 We already dis- 
cussed the presence of four peaks. The one nearest to the 
Fermi surface at negative energies is the one referred to 
by experimentalists when they refer to this renormaliza- 
tion. This band has a dispersion of order J (not shown 
on Fig. 26, but see Ref. 113). This result also agrees 
with VCPT as shown in Fig. 26 and CDMFT (Fig. 25a). 
As shown in Figs. 25a and 25d, whether the state is 
ordered or not the band width is similar. Analytical 
strong-coupling expansions 114 ' 115 and exact diagonaliza- 
tions also find the same result. To find detailed agree- 




FIG. 38: MDC at the Fermi energy for 10% hole-doped 
Ca 2 _ ;c Na ;c Cu02Cl2 from Ref. 85. 



ment with experiment, one needs to include t! and i". 116 
The evolution of the position of chemical potential for 
extremely small dopings as discussed in Ref. 117 is not 
reproduced by the strong-coupling calculations, although 
the result on chemical potential is somewhat material 
dependent. 118 

(ii) In the overdoped case, the Fermi surface is well 
defined. Although we have not shown any figures con- 
cerning this point, VCPT and CDMFT show the same 
result. 

(Hi) The evolution with doping of the electronic struc- 
ture has been mapped. It has shown the importance of 
antiferromagnetic correlations in the p-type underdoped 
cuprates and especially in the n-type ones in which the 
hot-spot physics is still observed at optimal doping. We 
will come back on the latter point for electron-doped 
cuprates in the following subsection. The strong-coupling 
results obtained with VCPT and CDMFT have a resolu- 
tion of order O.lt, which translates into about 30 meV. 
This is not enough to accurately measure the Fermi ve- 
locity, which was found to be doping independent in 
LSCO. 119 However, this suffices to compare with MDC 
curves obtained experimentally by integrating over an en- 
ergy range of about 60 meV, as shown in Fig. 38 obtained 
in Ref. 85 on Calcium oxyclorate Ca2- 2; Na :c Cu02Cl2, a 
10% hole-doped high temperature superconductor. The 
similarities between that figure and the CPT (Fig. 21), 
VCPT (Fig. 23) and CDMFT (Fig. 24) results is strik- 
ing. The agreement is better when no antiferromagnetic 
long-range order is assumed, as in the CPT case. The 
flattening of the band structure near {%, 0) observed ex- 
perimentally, can also be seen in CPT by comparing the 
top and middle EDC's taken at small and large U re- 
spectively on the left panel of Fig. 20. This flattening is 
associated with the pseudogap phenomenon. Recall that 
the theoretical results were obtained with t' = — 0.3t and 
t" = Q.2t. This in turn implies an electron-hole asym- 
metry that is observed experimentally. We come back to 
this in the following subsection. 

(iv) The overall d-wave symmetry of the superconduct- 
ing gap has been observed for both hole and electron dop- 
ing, supporting the universality of the pairing nature in 



the cuprates. In the next to next subsection, we discuss 
the phase diagram for competing antiferromagnetism and 
d-wave superconductivity and show striking similarities 
with the observations. 

(v) A normal-state pseudogap has been observed to 
open up at a temperature T* > T c in the underdoped 
regime with a d-wave form similar to the one of the su- 
perconducting gap. That statement is correct only in the 
hole-doped compounds. In electron-doped systems the 
pseudogap has a form that is not of d-wave shape. If 
T c comes from a universal pairing mechanism, a univer- 
sal mechanism may also be behind the pseudogap. As 
we have already discussed however, there are quantita- 
tive differences between strong and weak coupling mecha- 
nisms for both T c and the pseudogap. For electron-doped 
systems, we made quantitative predictions for the value 
of T* that have later been confirmed experimentally. All 
this is discussed further below. To date, in cluster meth- 
ods the pseudogap temperature has been studied only 
with DCA. 78 

(vi) A coherent quasiparticle peak below T c has been 
observed near (it, 0) whose spectral weight scales with the 
doping level x in the underdoped regime. We expect 
that it is a general result that long-range order will re- 
store quasiparticle like excitations in strongly correlated 
systems because gaps remove scattering channels near 
the Fermi level. This is clearly illustrated by comparing 
the upper and lower panels in Fig. 25 that contrast the 
same spectra with and without antiferromagnetic long- 
range order. We have not performed the analysis of our 
results yet that could tell us whether the spectral weight 
of the quasiparticle scales with x in the hole-underdoped 
regime. Our resolution may not be good enough to see 
the quasiparticle peak. Sharpening of the quasiparti- 
cle excitations in the superconducting state has however 
been observed in DCA. 120 

(vii) The presence of an energy scale of about 40 — 
80 meV in the quasiparticle dynamics manifests itself 
through a sharp dispersion renormalization and drop in 
the scattering rate observed at those energies at differ- 
ent momenta. In hole-doped systems there is a kink in 
the nodal direction that is already seen above T c while 
in the antinodal direction it appears only below T c . The 
energy scales and doping dependences of these two kinks 
are also different. 121 The energy resolution in VCPT and 
CDMFT is not sufficient to distinguish these subtleties. 
In electron-doped cuprates experiments 122 suggest that 
there is no observable kink feature, in agreement with 
the results presented in the following subsection. 



B. The pseudogap in electron-doped cuprates 

The ARPES spectrum of electron-doped cuprates is 
strikingly different from that of their hole-doped coun- 
terpart. The Fermi energy MDC's for the first quadrant 
of the Brillouin zone 86 are shown at the top of Fig. 39 for 
three different dopings. There is a very clear evolution 
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FIG. 39: Doping dependence of the MDC from experiments 
on NCCO with the corresponding EDC. From Ref. 86 



with doping. At the lowest dopings, there is no weight 
near (tt/2,tt/2), contrary to the hole-doped case shown 
in Fig. 38. For all dopings there is weight near (it, 0) in- 
stead of the pseudogap that appeared there in the hole- 
doped case. The EDC's, also shown on the bottom of 
Fig. 39, are drawn for a trajectory in the Brillouin zone 
that follows what would be the Fermi surface in the non- 
interacting case. Regions that are more green than red 
on the corresponding MDC's along that trajectory are re- 
ferred to as hot spots. On the EDC's we clearly see that 
hot spots do not correspond to simply a decrease in the 
quasiparticle weight Z. They truly originate from a pseu- 
dogap, in other words from the fact that the maximum 
is pushed away from zero energy. Even though the mea- 
surements are done at low temperature (T = 10 — 20K ) 
the energy resolution of about 60 meV makes the su- 
perconducting gap invisible. What is observed at this 
resolution is the pseudogap. 

The contrast between the location of the pseudogap in 
the hole and electron-doped compounds is clearly seen 
in Fig. 40 obtained from VCPT. 19 In that figure, the 
magnitude of the spectral weight is represented by the 
different colors as a function of frequency (in units of t) 
along different cuts of the Brillouin zone. In the bottom 
panel, for the hole-doped case, one observes the pseudo- 
gap near (w, 0). In the top panel, for the electron-doped 
case, it is only by zooming (inset) on the region for the 
Fermi energy crossing near (7r, 0) that one sees the d- 
wave superconducting gap. At 10% electron-doping, the 
pseudogap near (7r/2, tt/2) is apparent. In this case there 
is antiferromagnetic long-range order, but even if we use 
CPT that does not exhibit long-range order, there ap- 
pears a pseudogap in that region. 160 The main differ- 
ence between CPT and VCPT results is the bending back 
of the bands (for example around the symmetry axis at 
(ir/2, 7r/2)) caused by halving of the size of the Brillouin 
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FIG. 40: Intensity plot of the spectral function as a function 
of lu in units of t and wave vector from VCPT for U = 8t 
t' = -O.St, t" = 0.2t and n = 0.93 at the bottom and n = 1.10 
(electron-doped) at the top. The Lorentzian broadening is 
0.12t in the main figure and 0.04t in the inset that displays 
the d-wave gap. Top panel is for the electron-doped case in 
the right-hand panel of Fig. 23, while bottom panel is for the 
hole-doped case on the left of Fig. 23. From Ref. 19. 
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FIG. 41: Experimental Fermi surface plot (MDC at the Fermi 
energy) for NCCO (left) and corresponding energy distribu- 
tion curves (right) for 15% electron-doping. From Ref. 86 

zone in the antiferromagnetic case. Form factors are 
such that the intensity is not symmetric even if the dis- 
persion is. The faint band located at an energy about t 
below the Fermi energy near (7r/2, tt/2) was also found in 
Ref. 89 by a one-loop spin-wave calculation around the 
Hartree-Fock antiferromagnetic ground state at U = 8t. 
Experimentalists 86 have suggested the existence of these 
states. The VCPT results go well beyond the spin-wave 
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FIG. 42: EDC A<(k,o;) = A(k,u)f(u>) along the Fermi sur- 
face calculated in TPSC (left) at optimal doping for t' = 
-0.175i, t" = 0.05t, t = 350 meV and corresponding ARPES 
data on NCCO (right). From Ref. 87. 




FIG. 43: EDC A < (k,uj) = A(k,w)/(w) along two other di- 
rections calculated for t' = -0.175t, t" = 0.05i, t = 350 meV 
in TPSC (left column) and corresponding ARPES data on 
NCCO (right column). From Ref. 87. 



calculation (dashed lines in Fig. 40) since one can also 
see numerous features in addition to remnants of the lo- 
calized atomic levels around +5t and — lOt. 

The optimally doped case is the real challenge for 
strong-coupling calculations. The spin-wave approach in 
Ref. 89 never shows the weight near (ir/2,ir/2) that is 
seen in experiment (Fig. 41). Early mean-field calcula- 
tions by Kusko et al. 90 suggest that this (ir/2,ir/2) fea- 
ture appears for U — it. This is very small compared 
with U of the order of the bandwidth 8t necessary to 
have a Mott insulator at half-filling. We already dis- 
cussed in Sec. Ill A that both CPT and TPSC show 
that weight near (tt/2,it/2) appears for U not too large, 
say of order 6t. This same result is also obtained in the 
Kotliar-Ruckenstein slave boson approach. 88 

Since TPSC is valid for a system of infinite size, 
we present detailed comparisons 87 with experiment 86 
on Ndi.ssCeo.isCuC^, an electron-doped cuprate. We 
take t' = -0.175i, t" = 0.05i. Results obtained with 
t' = —0.275 are very close to those we present. With the 
values used in CPT, t' = -0M, t' = 0.2t, U = 6i, TPSC 
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FIG. 44: EDC 4<(k,w) = A(k,w)f(ui) along the Fermi sur- 
face shown in the insets for (b) n = 1.10, U = 6.25t. Lines 
are shifted by a constant for clarity. From Ref. 87. 




does not lead to strong enough antiferromagnetic fluc- 
tuations to obtain non-trivial effects in the temperature 
range studied, (3 = 20t. We take t = 350 meV. Fig. 39 
shows the correspondence between EDC and MDC. Com- 
parisons with experimental EDC at wave vectors along 
the non-interacting Fermi surface appear in Fig. 42 for 
U = 5.75i and 15% doping (n = 1.15). The dashed lines 
indicate the quite detailed agreement between theory and 
experiment. At the hot spot, (middle dashed line), the 
weight is pushed back about 0.2 eV and there is a very 
small peak left at the Fermi surface, as in the experiment. 
If U is not large enough the antiferromagnetic fluctua- 
tions are not strong enough to lead to a pseudogap. As in 
CPT (Fig. 21), if U is too large the (tt/2, tt/2) weight dis- 
appears, as illustrated earlier in Fig. 29. In Fig. 43, cuts 
along the (0,0) to {n,w) and (0.657r,0) to (0.657T, tt) di- 
rections are compared with experiment. Again the peak 
positions and widths are very close, except for some ex- 
perimental tails extending in the large binding energy 
direction. The theoretical results have similar asymme- 
try, but not as pronounced. Experimentally, the large 
binding energy tails ("the background") are the least 
reproducible features from sample to sample, especially 
for wave vectors near the Fermi surface". 161 The experi- 
mental renormalized Fermi velocities are 3.31 x 10 5 m/s 
and 3.09 x 10 5 m/s along the zone diagonal and along 
the (n, 0)-(7T,7r) direction, respectively. The correspond- 
ing renormalized Fermi velocities obtained by TPSC are 
3.27 x 10 5 m/s and 2.49 x 10 5 m/s, respectively. The 
agreement is very good, particularly along the diagonal 
direction. The bare Fermi velocities are renormalized in 
TPSC by roughly a factor of two. 123 

As we move towards half-filling, we have to increase U 
slightly to find agreement with experiment, as discussed 
earlier in Fig. 29. Fig. 44 shows how well the EDC's agree 
for a Fermi surface cut at 10% doping (n = 1.10). The 
increase is expected physically from the fact that with 
fewer electrons the contribution to screening that comes 
from Thomas Fermi physics should not be as good. This 
is also consistent with the fact that a larger value of U 
is necessary to explain the Mott insulator at half-filling. 
It would also be possible to mimic the ARPES spectrum 
by keeping U fixed and changing the hopping parameters, 



FIG. 45: Hot spots from quasi-static scatterings off antiferro- 
magnetic fluctuations (renormalized classical regime). 
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FIG. 46: Semi-log plot of the AFM correlation length (in units 
of the lattice constant) against inverse temperature (in units 
of J = 125 meV). Filled symbols denote calculated results 
and empty ones experimental data of Ref. 124 and Ref. 125 
(ac = 0.15). From Ref. 87. 



but the changes would be of order 20%, which does not 
appear realistic. 87 

We have already explained that the physics behind the 
pseudogap in TPSC is scattering by nearly critical anti- 
ferromagnetic fluctuations. This is illustrated in Fig. 45. 
If this explanation is correct, the antiferromagnetic cor- 
relation length measured by neutron scattering should 
be quite large. The results of the measurement 124,125 
and of the TPSC calculations are shown in Fig. 46. The 
agreement is again surprisingly good. As we move to 
smaller dopings n = 1.1 (not shown) the agreement be- 
comes less good, but we do expect TPSC to deteriorate 
as U increases with underdoping. The arrow points to 
the temperature where EDC's shown earlier were calcu- 
lated. Note however that the neutron measurements were 
done on samples that were not reduced, by contrast with 
the ARPES measurements mentioned earlier. We are ex- 
pecting experiments on this subject. 162 . We should point 
out that the EDC's depend strongly on temperature and 
on the actual value of the antiferromagnetic correlation 
length only in the vicinity of the temperature where there 
is a crossover to the pseudogap regime. Decreasing the 
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FIG. 47: Pseudogap temperature T* (filled circles de- 
note T* calculated from TPSC, empty ones experimental 
data extracted from optical conductivity. 127 ) Empty trian- 
gles are experimental Neel temperatures Tn ■ The samples 
are reduced. 124 From Ref. 87. 



FIG. 48: The generic phase diagram of high-T c superconduc- 
tors, from Ref. 110. There should also be a pseudogap line on 
the electron-doped side. It was not well studied at the time 
of publication of that paper. 



temperature makes the w = peaks near (ir, 0) sharper 87 
as observed experimentally. 126 

The ARPES pseudogap temperature T* has been pre- 
dicted with TPSC. 87 The predictions are shown by the 
solid line in Fig. 47. The pseudogap temperature ob- 
served in optical experiments 127 is shown by the open 
circles. It differs from the ARPES result, especially as 
we move towards optimal doping. The size of the pseu- 
dogap observed in the optical experiments 127 (10T*) is 
comparable to the ARPES pseudogap. The solid line in 
Fig. 47 contains several predictions. If we look at 13% 
doping (n = 1.13), the line predicts T* ~ 250JT. Ex- 
periments that were done without being aware of this 
prediction 128 have verified it. It would be most inter- 
esting to do neutron scattering experiments on the same 
samples to check whether the antifcrromagnetic corre- 
lation length £ and the thermal de Broglie wave length 
tlth are comparable at that temperature, as predicted by 
TPSC. Fig. 47 also predicts that the pseudogap induced 
by antifcrromagnetic fluctuations will disappear at the 
quantum critical point where long-range antiferromag- 
netic order disappears, in other words it will coincide 
with the crossing of the experimentally observed Neel 
temperature (dashed line with triangles in Fig. 47) with 
the zero temperature axis (if that crossing is not masked 
by the superconducting transition). Recent optical con- 
ductivity experiments 129,130 confirm this prediction as 
well. 

In TPSC, superconducting fluctuations can also lead 
to a pseudogap by an analogous mechanism. 23 
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FIG. 49: Antiferromagnetic order parameter m (dashed) and 
d-wave (solid) order parameter obtained from CDMFT on a 
2x2 cluster. The result obtained by forcing m = is also 
shown as a thin dashed line. 



half-filling. There are other features on the phase di- 
agram, in particular checkerboard patterns 131 or stripe 
phases 132 that appear in general close to the region where 
antifcrromagnetism and superconductivity come close to 
each other. Before we try to understand these more de- 
tailed features, one should understand the most impor- 
tant phases. In the previous subsection we have discussed 
the pseudogap phase, in particular on the electron-doped 
side (not indicated on Fig. 48). A recent review of the 
pseudogap appears in Ref. 133. In the following, we dis- 
cuss in turn the phase diagram and then the nature of 
the superconducting phase itself and its relation to the 
Mott phenomenon. 



C. The phase diagram for high-temperature 
superconductors 

The main features appearing in the phase diagram 
of high-temperature superconductors are the pseudogap 
phase, the antiferromagnetic phase and the d-wave su- 
perconducting phase. Fig. 48 110 shows the typical dia- 
gram with hole doping to the right and electron doping 
to the left. Zero on the horizontal axis corresponds to 



1. Competition between antiferromagnetism and 
superconductivity 

We have already shown in Fig. 32 the prediction of 
VCPT for the zero-temperature phase diagram. 19 Here, 
we just point out how closely the position of the anti- 
ferromagnetic phase boundary, appearing in the lower 
panel, coincides with the experimental phase diagram in 



26 



Fig. 48. (Note that electron concentration increases from 
right to left on this experimental phase diagram) . In par- 
ticular, there is little size dependence to the position of 
this boundary, (6 to 10 sites) and in addition the depen- 
dence on the value of U is also weak, as can be seen from 
Fig. 34. Hence, the positions of the antifcrromagnctic 
phase boundaries is a robust prediction of VCPT. The 
CDMFT result for a four site cluster in a bath is shown 
in Fig. 49 for t' = -0.3t, t" = and U = 8t. The agree- 
ment with experiment is not as good. Despite the useful 
presence of a bath in CDMFT, the cluster itself is of size 
2x2, which is probably smaller than the Cooper pair 
size. We can obtain results closer to those of VCPT by 
increasing the variational space. 

The d-wave superconducting order parameter on the 
top panel of Fig. 32 shows more size dependence than 
the antiferromagnetic order parameter. Nevertheless, 
there are some clear tendencies: (a) d-wave supercon- 
ductivity can exist by itself, without antiferromagnetism. 
The vertical lines indicate the location of the end of 
the antiferromagnetic phase for the various system sizes 
to help this observation, (b) The range where d-wave- 
superconductivity exists without antiferromagnetism, is 
about three times larger on the hole than on the electron- 
doped side, as observed experimentally, (c) As system 
size increases, the maximum d-wave order parameter 
is larger on the hole than on the electron-doped side, 
(d) The tendency to have coexisting antiferromagnetism 
and d-wave superconductivity is rather strong on the 
electron-doped side of the phase diagram. This is ob- 
served experimentally 134 but only over a rather narrow 
region near optimal doping. Recent experiments 135 chal- 
lenge this result, others 136,137 indicate that antiferromag- 
netism can be induced from the d-wave superconducting 
phases with very small fields, (e) On the hole-doped side, 
d-wave superconductivity and antiferromagnetism coex- 
ist for a very narrow range of dopings for system size 
N c = 6, for a broad range extending to half- filling for 
N c = 8 and not at all for N c = 10. In other words, 
the tendency to coexistence is not even monotonic. We 
interpret this result as a reflection of the tendency to 
form stripes observed experimentally on the hole-doped 
side. 132 ' 138 ' 139 We cannot study systems large enough 
to allow for striped inhomogeneous states to check this 
statement. 

The more realistic two-band model has also been stud- 
ied using DCA. 104 The results are shown on Fig. 50. 
Electron concentration increases from right to left. This 
phase diagram is very close to that obtained with the 
same method from the one-band Hubbard model 104 with 
t' = -0.3t, t" = 0, U = St. The qualitative results agree 
with the other calculations and with experiment: anti- 
ferromagnetism extends over a narrower doping range for 
hole than for electron doping and d-wave superconductiv- 
ity by itself exists over a broader range for the hole-doped 
case than for the electron-doped case. The actual ranges 
where antiferromagnetism and d-wave superconductivity 
exist are not in as good an agreement with experiment 




FIG. 50: Phase diagram obtained from DCA for U = 8t for 
the two-band model. From Ref. 104 



as in the VCPT case. However, as in CDMFT, the sys- 
tem sizes, 2x2, arc very small. Overall then, quantum 
cluster methods, VCPT in particular, allow us to obtain 
from the Hubbard model the two main phases, antifer- 
romagnetic and d-wave superconducting, essentially in 
the observed doping range of the zero-temperature phase 
diagram. At finite temperature, DCA and TPSC agree 
on the value of T c for the particle-hole symmetric model 
at 10% doping and U = At. Recent studies of the irre- 
ducible vertex using DCA 140 also show that in the weak- 
coupling limit the particle-particle d-wave channel leads 
to an instability driven by antiferromagnetic fluctuations 
as temperature decreases, as found in TPSC. 

To understand the effect of pressure on the phase dia- 
gram, note that U/t should decrease as pressure increases 
since the increase in the overlap between orbitals should 
lead mainly to an increase in t. Hence, as can be de- 
duced from Fig. 37, applying pressure should lead to a 
decrease in the value of T c at weak coupling, concomi- 
tant with the decrease in antiferromagnetic fluctuations 
that lead to pairing in the weak coupling case. This is 
indeed what pressure does experimentally in the case of 
electron-doped high-temperature superconductors, 141 re- 
inforcing our argument that near optimal doping they are 
more weakly coupled. It is widely known on the other 
hand that pressure increases T c in hole-doped systems. 
That is consistent with the strong-coupling result that 
we found in VCPT and CDMFT, namely that the max- 
imum d-wave order parameter in that case scales with 
J = 4t 2 /U, a quantity that increases with t and hence 
pressure. Whereas in the weak coupling case supercon- 
ductivity is a secondary phenomenon that occurs after 
antiferromagnetic fluctuations have built up, in strong 
coupling they can be two distinct phenomena as can be 
seen from the phase diagram, even though they arise from 
the same microscopic exchange interaction represented 
by J. 
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2. Anomalous superconductivity near the Mott transition 

Superconductivity in the undcrdoped regime is very 
much non-BCS. First of all, we notice in Fig. 35 ob- 
tained in CDMFT 17 that at strong coupling the d-wave 
superconducting order parameter vanishes as we move 
towards half-filling even in the absence of long-range an- 
tiferromagnetic order. In other words, the Mott phe- 
nomenon by itself suffices to destroy d-wave superconduc- 
tivity. This conclusion is reinforced by the fact that at 
weak coupling (U = At) where there is no Mott localiza- 
tion, d-wave superconductivity survives at half-filling. In 
the presence of antifcrromagnetic long-range order, that 
last statement would not be true, as confirmed by VCPT 
calculations in Fig. 36: at U — At d-wave superconduc- 
tivity survives at half-filling if we do not allow for an- 
tifcrromagnetic long-range order but it disappears if we 
do. In BCS theory, the presence of an interaction J that 
leads to attraction in the d-wave channel would lead at 
T = to d-wave superconductivity at all dopings includ- 
ing half-filling, unless we allow for competing long-range 
order. At strong coupling, no long-range order is neces- 
sary to destroy d-wave superconductivity. 

Superconductivity at strong coupling 142,143 also differs 
from BCS in the origin of the condensation energy. Sup- 
pose we do BCS theory on the attractive Hubbard model. 
Then, as in the usual BCS model, kinetic energy is in- 
creased in the superconducting state because the Fermi 
surface is no-longer sharp. On the other hand, in the 
superconducting phase there is a gain in potential en- 
ergy. The reverse is true at strong coupling. This re- 
sult follows from DCA 120 and is in agreement with the 
kinetic energy drop in the superconducting state that 
has been estimated from the /-sum rule in optical con- 
ductivity experiments. 144-146 Photoemission data 147 had 
also suggested this kinetic energy drop in the supercon- 
ducting state. A crossover from non-BCS-like to BCS 
behavior in the condensation mechanism as we go from 
underdoping to overdoping has also been seen recently 
experimentally. 144 We do not seem to have the resolution 
to find that crossover since the condensation energy be- 
comes very small on the overdoped side. We expect that 
crossover from strong to weak coupling will also lead to a 
change from a kinetic-energy driven to a potential-energy 
driven pairing mechanism. This is confirmed by CDMFT 
calculations for the attractive Hubbard model. 148 

A third way in which superconductivity in the under- 
doped regime is non-BCS is that the drop in the order 
parameter as we go towards half-filling is accompanied by 
an increase in the gap as measured in the single-particle 
density of states. Fig. 6 of Ref. 149 summarizes the ex- 
perimental evidence for the increase in the size of the gap. 
That increase, observed in the CDMFT calculation of the 
gap, is illustrated in Fig. 51. 17 That gap has essentially 
the same size as that observed in the normal pseudogap 
state 67 . 
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FIG. 51: The gap in the density of states of the dSC as a 
function of filling for U = 8t, t' = —0.3t as calculated in 
CDMFT on a 2 x 2 cluster. From Ref. 17. 



V. CONCLUSION, OPEN PROBLEMS 

High-temperature superconductivity has forced both 
experimentalists and theorists to refine their tools and to 
develop new ones to solve the puzzles offered by this re- 
markable phenomenon. From a theoretical perspective, 
the original suggestion of Anderson 1 that the physics 
was in the one-band Hubbard model is being confirmed. 
In the absence of ab initio methods to tell us what is 
the correct starting point, such insight is essential. The 
non-perturbative nature of the phenomenon has however 
forced theorists to be extremely critical of each other's 
theories since none of them can pretend that a small pa- 
rameter controls the accuracy of the approximations. 

If theorists are to convince each other and experimen- 
talists that a solution of the high-temperature super- 
conductivity problem has been found, then the theories 
have to give quantitative results and to make predic- 
tions. Unlike most traditional problems in condensed 
matter physics however, the non-perturbative nature of 
the problem means that no simple mean-field like the- 
ory can be trusted, even if it seems to agree qualitatively 
with experiment. In fact several such theories have been 
proposed 13,14 ' 150 ' 151 not long after the experimental dis- 
covery of the phenomenon but they have not been ac- 
cepted immediately. Theories have to be internally con- 
sistent, they have to agree with exact results whenever 
they are available, and then they can be compared with 
experiments. If there is a disagreement with experiment, 
the starting point (one-band Hubbard model) needs to 
be reconsidered. When approaches developed on the 
basis of weak-coupling ideas agree at intermediate cou- 
pling with approaches developed on the basis of strong- 
coupling ideas, then one gains confidence in the validity 
of the results. We have argued that such concordance is 
now found in a number of cases and that corresponding 
rather detailed quantitative agreement with experiment 
can be found. In a non-perturbative context it becomes 
essential to also cross check various approaches. 

The main theoretical methods that we have dis- 
cussed are those that we have developed or perfected 
or simply used in our group: The Two-Particle Self- 
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Consistent approach that is based on weak-coupling but 
non-perturbative ideas (no diagrams are involved), as 
well as heavily numerical approaches such as QMC and 
various quantum cluster methods, VCPT and CDMFT. 

Based on our own work and that of many others, 
we think the following experimental facts about high- 
temperature superconductivity can be reproduced very 
accurately by calculations for the one-band Hubbard 
model with U in the intermediate coupling range (U ~ 
8f) with t ~ 350 mcV, and hopping parameters t' 
and t" close to the values suggested by band structure 
calculations, 103 namely t' = -0.3t, t" = 0.2t. 

(i) In the one-band Hubbard model the main phases of 
the zero-temperature phase diagram, namely antiferro- 
magnetic and d-wave superconducting, appear very near 
the observed ranges for both the hole- and electron-doped 
cases. 

(ii) The normal state is unstable to a d-wave supercon- 
ducting phase in a temperature range that has the correct 
order of magnitude. As usual the value of T c is the most 
difficult quantity to evaluate since one must take into ac- 
count Kostcrlitz-Thouless physics as well as the effect of 
higher dimensions etc, so this level of agreement must be 
considered satisfying. 

(iii) The ARPES MDC at the Fermi energy and 
the EDC near the Fermi energy are qualitatively well 
explained by cluster calculations for both hole- and 
electron-doped cases. These comparisons, made at a res- 
olution of about 30 to 60 meV are not very sensitive to 
long-range order, although order does influence the re- 
sults. One is mainly sensitive to the pseudogap, so this 
is the main phenomenon that comes out from the model. 
Energy resolution is not good enough to see a kink. More 
details about what aspects of ARPES arc understood 
may be found in Sec. IV A. 

(iv) In the case of electron-doped cuprates, the value of 
U near optimal doping seems to be in the range U ~ 6i, 
which means that it is accessible to studies with TPSC 
that have better resolution. In that case, the agreement 
with experiment is very accurate, even if there is room 
for improvement and a need for further experiments. In 
addition, the value of T* for 13% doping has been pre- 
dicted theoretically before it was observed experimen- 
tally, one of the very rare predictions in the field of 
high-temperature superconductors. All of this agreement 
with ARPES data is strong indication that U <~ 6t is 
appropriate to describe electron-doped superconductors 
near optimal doping. Additional arguments come from 
the pressure dependence of the superconducting transi- 
tion temperature T c , which increases with t/U contrary 
to the strong-coupling result, and from simple ideas on 
Thomas-Fermi screening. The latter would predict that 
the screened interaction scales like dfi/dn and CPT re- 
sults do lead to d/i/dn smaller on the electron- than on 
the hole-doped side. 77 In addition, the optical gap at half- 
filling is smaller in electron- than in hole-doped systems. 

What is the physics? The physics of the antiferro- 
magnetic phase at both weak and strong coupling is 



well understood and needs no further comment. For the 
pseudogap, we have argued that there seems to be two 
mechanisms, a weak coupling one that involves scattering 
off critical fluctuations and that is very well understood 
within TPSC, and a strong-coupling one where there is 
no need for large correlation lengths. There is no sim- 
ple physical picture for the latter mechanism although 
the fact that it does not scale with J but with t seems 
to suggest forbidden hopping. The pseudogap is clearly 
different from the Mott gap. Whether there is a phase 
transition as a function of U that separates the weak 
and strong coupling regimes or whether there is only a 
crossover is an open question. The shape of the MDC's 
at the Fermi energy clearly show in any case that in some 
directions wave vector is not such a bad quantum num- 
ber whereas in the pseudogap direction, a "localized" or 
"almost localized" particle-like picture would be appro- 
priate. In fact the pseudogap occurs near the intersection 
with the antiferromagnetic zone boundary that turns out 
to also be the place where umklapp processes are possi- 
ble. In other words, the presence of a lattice is extremely 
important for the appearance of the pseudogap. We have 
seen that with spherical Fermi surfaces the Fermi liquid 
survives even for large U. The dichotomy between the 
wave description inherent to the Fermi liquid and the 
particle (localized) description inherent to the Mott phe- 
nomenon seems to be resolved in the pseudogap phase by 
having certain directions where electrons are more wave- 
like and other directions where particle-like (gapped) be- 
havior appears. The latter behavior appears near regions 
where the presence of the lattice is felt through umklapp 
processes. 

It is clear that when weak-coupling-like ideas of quasi- 
particles scattering off each other and off collective ex- 
citations do not apply, a simple physical description be- 
comes difficult. In fact, knowing the exact wave functions 
would give us the solution but we would not know how 
to understand "physically" the results. 

This lack of simple physical images and the neces- 
sity to develop a new discourse is quite apparent for 
d-wave superconductivity. At weak coupling exchange 
of slow antiferromagnetic fluctuations is at the origin of 
the phenomenon, while at strong-coupling the fact that 
the maximum value of the d-wave order parameter scales 
with J tells us that this microscopic coupling is impor- 
tant, even though there is no apparent boson exchange. 
This is where mean-field like theories 13,14 ' 150 or varia- 
tional approaches 4 can help when they turn out to give 
results that are confirmed by more accurate and less bi- 
ased methods. 

There are many open problems, some of which are ma- 
terial dependent and hence may depend on interactions 
not included in the simplest Hubbard model. We have 
already mentioned the problem of the chemical potential 
shift in ARPES for very small dopings 117 that seems to 
be somewhat material dependent 118 . It would also be im- 
portant to understand additional inhomogeneous phases 
that are observed in certain high-temperature supercon- 



29 



ductors. That is extremely challenging for quantum clus- 
ter methods and unlikely to be possible in the very near 
future, except for inhomogencitics of very short wave 
length. Also, we still need to improve concordance be- 
tween the methods before we can make predictions that 
are quantitative at the few percent level for all physical 
quantities. Apart from DCA, there are no quantum clus- 
ter methods that have been developed yet to study two- 
particle response functions that are necessary to obtain 
results on the superfluid density and on transport in gen- 
eral. Transport studies are being completed in TPSC. 152 
Such studies are crucial since they are needed to answer 
questions such as: (i) Why is it that for transport proper- 
ties, such as optical conductivity, the number of carriers 
appears to scale with doping whereas in ARPES the sur- 
face of the Brillouin zone enclosed by the apparent Fermi 
surface appears to scale with the number of electrons? 
Is it because the weight of quasiparticles at the Fermi 
surface scales like the doping or because of vertex correc- 
tions or because of both? (ii) Can we explain a vanishing 
superfluid density as doping goes to zero 153 only through 
Mott physics or can competing order do the job. 163 , 94 
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APPENDIX A: LIST OF ACRONYMS 



After twenty years all the problems are not solved, but 
we think that we can say with confidence that the es- 
sential physics of the problem of high-temperature su- 
perconductivity is in the one-band Hubbard model. At 
least the pseudogap, the antiferromagnetic and the d- 
wave superconducting phases come out from the model. 
Refinements of that model may however be necessary as 
we understand more and more details of the material- 
specific experimental results. 

Has a revolution been necessary to understand the ba- 
sic physics of high-temperature superconductors? Cer- 
tainly, it has been necessary to change our attitude to- 
wards methods of solution. We have seen that to study 
intermediate coupling, even starting from weak coupling, 
it has been necessary to drop diagrams and to rely in- 
stead on sum rules and other exact results to devise a 
non-perturbative approach. At strong coupling we had to 
accept that numerical methods are essential for progress 
and that we need to abandon some of the traditional 
physical explanations of the phenomena in terms of el- 
ementary excitations. Even though progress has been 
relatively slow, the pace is accelerating in the last few 
years and there is hope that in a few years the problem 
will be considered for the most part solved. The the- 
oretical methods (numerical and analytical) that have 
been developed and that still need to be developed will 
likely remain in the tool box of the theoretical physicist 
and will probably be useful to understand and perhaps 
even design other yet undiscovered materials with inter- 
esting properties. The success will have been the result 
of the patient and focused effort of a large community of 
scientists fascinated by the remarkable phenomenon of 
high- temper at ure superconductivity. 



ARPES: Angle Resolved Photoemission Spec- 
troscopy: Experiment from which one can 
extract A(k,w) / (w) . 

CPT: Cluster Perturbation Theory: Cluster 
method based on strong coupling perturbation 
theory. 52 ' 53 - 62 

CDMFT: Cellular Dynamical Mean Field Theory: A 
cluster generalization of DMFT that allows one to 
take into account both wave vector and frequency 
dependence of the self-energy. 58 It is best formu- 
lated in real space. 

DCA: Dynamical Cluster approximation: A cluster gen- 
eralization of DMFT that allows one to take into ac- 
count both wave vector and frequency dependence 
of the self-energy based on coarse graining of the 
self-energy in reciprocal space. 55,154 

DMFT: Dynamical Mean Field Theory: This approach 
is exact in infinite dimension. It takes the frequency 
dependence of the self-energy into account and in- 
cludes both the Mott and the Fermi liquid limits. 2,3 

EDC: Energy Dispersion Curves: A representation of 
A (k,w) / (w) at fixed k as a function of uj. 

FLEX: Fluctuation Exchange Approximation: A con- 
serving many-body approach, similar in spirit to 
Eliashbcrg theory. 22 

MDC: Momentum Dispersion Curves: A representation 
of A(k,u>)f((jj) at fixed lo as a function of k. 

QMC: Quantum Monte Carlo: Determinental 
approach 155 . This provides an essentially ex- 
act solution to the model for a given system size 
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and within statistical errors that can be made 
smaller by performing more measurements. 

RPA: Random Phase Approximation. 

TPSC: Two-Particle Self-Consistent Approach: Based 
on sum rules and other constraints, allows to treat 
the Hubbard model non-perturbatively in the weak 
to intermediate coupling limit. 23 ' 34 



VCA: Variational Cluster Approach. Analogous to 
CDMFT but with a convergence criterion based on 
an extremum principle. In the applications quoted 
here, there is no bath. 56 - 64 - 73 - 156 

VCPT: Variational Cluster Perturbation Theory. In the 
present paper synonymous with VCA. 



1 P. Anderson, Science 235, 1196 (1987). 

2 A. Georges, G. Kotliar, W. Krauth and M. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

3 A. Georges and G. Kotliar, Phys. Rev. B 45, 6479 (1992); 
M. Jarrell, Phys. Rev. Lett. 69, 168 (1992). 

4 A. Paramekanti, M. Randeria and N. Trivedi, Phys. 
Rev. B 70, 54504 (2004). 

5 J. Hirsch and H. Lin, Phys. Rev. B 37, 5070 (1988). 

6 M. Beal-Monod, C. Bourbonnais and V. Emery, Phys. 
Rev. B 34, 7716 (1986). 

7 K. Miyake, S. Schmitt-Rink and C. Varma, Phys. Rev. 
B 34, 6554 (1986). 

8 D. Scalapino, E. Loh and J. Hirsch, Phys. Rev. B 34, 
8190 (1986). 

9 W. Kohn and J. M. Luttinger, Phys. Rev. Lett. 15, 524 
(1965). 

10 A. Leggett, Rev. Mod. Phys. 47, 331 (1975). 
45 N. Bulut, D.J. Scalapino, SR. White, Phys. Rev. B 47, 
14 599 (1993). 

12 S. Zhang, J. Carlson and J. E. Gubernatis, Phys. Rev. 
Lett. 78, 4486 (1997). 

13 M. Inui, S. Doniach, P. Hirschfeld and A. Ruckenstein, 
Phys. Rev. B 37, 2320 (1988). 

14 G. Kotliar and J. Liu, Phys. Rev. Lett. 61, 1784 (1988). 

15 T. Maier, M. Jarrell, T. Pruschke and J. Keller, Phys. 
Rev. Lett. 85, 1524 (2000). 

16 D. Poilblanc and D. Scalapino, Phys. Rev. B 66, 052513 
(2002). 

17 S. S. Kancharla, M. Civelli, M. Capone, B. Kyung, 
D. Senechal, G. Kotliar and A.-M. Tremblay, Anomalous 
superconductivity in doped Mott insulators (2005), cond- 
mat/0508205. 

18 T. Maier, M. Jarrell, T. Schulthess, P. Kent and 
J. White, hys. Rev. Lett. 95, 237001 (2005). 

19 D. Senechal, P.-L. Lavertu, M.-A. Marois and A.-M. S. 
Tremblay, Phys. Rev. Lett. 94, 156404 (2005). 

20 T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 
netism (Springer-Verlag, 1985). 

21 R. S. Markiewicz, Phys. Rev. B 70, 174518 (2004). 

22 N. Bickers and D. Scalapino, Ann. Phys. (USA)193, 206 
(1989). 

23 Y. Vilk and A.-M. Tremblay, J. Phys I (France) 7, 1309 
(1997). 

24 Y. Vilk and A.-M. Tremblay, Europhys. Lett. 33, 159 
(1996). 

25 N. Furukawa, T. M. Rice and M. Salmhofer, Phys. Rev. 
Lett. 81, 3195 (1998). 

26 C. Halboth and W. Metzner, Phys. Rev. B 61, 7364 
(2000). 

27 V. Hankevych and F. Wegner, Eur. Phys. J. B 31, 333 



(2003). 

28 C. Honerkamp and M. Salmhofer, Phys. Rev. Lett. 87, 
187004 (2001). 

29 F. Wegner and V. Hankevych, Acta Phys. Pol. B 34, 
497 (2003). 

30 A. A. Katanin and A. P. Kampf, Phys. Rev. Lett. 93, 
106406 (2004). 

31 D. Zanchi, Europhys. Lett. 55, 376 (2001). 

32 G. Mahan, Many-Particle Physics, 3rd edition, Section 
6.4.4 (Kluwer/Plenum, 2000). 

33 Y. Vilk and A.-M. Tremblay, J. Phys. Chem. Solids (UK) 
56, 1769 (1995). 

34 S. Allen, A.-M. Tremblay and Y. M. Vilk, in Theoreti- 
cal Methods for Strongly Correlated Electrons, edited by 
D. Senechal, C. Bourbonnais and A.-M. Tremblay, CRM 
Series in Mathematical Physics, Springer (2003). 

35 G. Baym, Physical review 127, 1391 (1962). 

36 S. Moukouri, S. Allen, F. Lemay, B. Kyung, D. Poulin, 
Y. Vilk and A.-M. Tremblay, Phys. Rev. B 61, 7887 

(2000) . 

37 Y. Vilk, L. Chen and A.-M. Tremblay, Phys. Rev. B 49, 
13267 (1994). 

38 S. Roy, C. Brillon and A.-M. Tremblay, unpublished. 

39 S. Allen, B. Kyung and A.-M. Tremblay, unpublished. 

40 B. Kyung, J.-S. Landry and A. M. S. Tremblay, Phys. 
Rev. B 68, 174502 (2003). 

41 S. Allen and A.-M. Tremblay, Phys. Rev. B 64, 075115 

(2001) . 

42 L. Lilly, A. Muramatsu and W. Hanke, Phys. Rev. Lett. 
65, 1379 (1990). 

43 B. Kyung, S. Allen and A.-M. Tremblay, Phys. Rev. B 
64, 075116 (2001). 

44 S. White, D. J. Scalapino and R. Sugar, Phys. Rev. B 
39, 839 (1989). 

45 N. Bulut, D. Scalapino and S. White, Phys. Rev. B 47, 
2742 (1993). 

46 A. Dare, Y. Vilk and A.-M. Tremblay, Phys. Rev. B 53, 
14236 (1996). 

47 L. Chen, C. Bourbonnais, T. Li and A.-M.-S. Tremblay, 
Phys. Rev. Lett. 66, 369 (1991). 

48 V. Hankevych, B. Kyung and A.-M. Tremblay, Phys. 
Rev. B 68, 214405 (2003). 

49 N. Trivedi and M. Randeria, Phys. Rev. Lett. 75, 312 
(1995). 

50 A. Moreo, D. J. Scalapino and S. R. White, Phys. Rev. 
B 45, 7544 (1992). 

51 G . Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. Parcollet, C.A. Marianetti, Electronic Structure Cal- 
culations with Dynamical Mean-Field Theory: A Spectral 
Density Functional Approach (2005), cond-mat/0511085. 



31 



52 C. Gros and R. Valenti, Phys. Rev. B 48, 418 (1993). 
j3 D. Senechal, D. Perez and M. Pioro-Ladriere, Phys. Rev. 
Lett. 84, 522 (2000). 

54 A. Schiller and K. Ingersent Phys. Rev. Lett. 75, 113 
(1975) 

55 M. H. Hettler, A. N. Tahvildar-Zadeh, M. Jarrell, T. Pr- 
uschke and H. Krishnamurthy, Phys. Rev. B 58, R7475 
(1998). 

56 M. Potthoff, Eur. Phys. J. B 36, 335 (2003). 

57 M. Potthoff, M. Aichhorn and C. Dahnken, Phys. Rev. 
Lett. 91, 206402 (2003). 

58 G. Kotliar, S. Savrasov, G. Palsson and G. Biroli, Phys. 
Rev. Lett. 87, 186401 (2001). 

59 T. Stanescu and G. Kotliar, Cluster Dynamical Mean 
Field Theories: A strong coupling perspective (2005), 
cond-mat/0508302. 

60 T. D. Stanescu and G. Kotliar, Phys. Rev. B 70, 205112 

(2004) . 

61 T. Maier, M. Jarrell, T. Pruschke and M. H. Hettler, 
Rev. Mod. Phys. 77, 1027 (2005). 

62 D. Senechal, D. Perez and D. Plouffe, Phys. Rev. B 66, 
075129 (2002). 

63 M. Potthoff, Non-perturbative construction of the 
Luttinger-Ward functional (2004), cond-mat/0406671. 

64 M. Potthoff, Adv. Solid State Phys. 45, 135 (2005). 

65 M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 
(1994). 

66 C. J. Bolech, S. S. Kancharla and G. Kotliar, Phys. Rev. 
B 67, 075110 (2003). 

67 B. Kyung, S. Kancharla, D. Senechal, A. M. S. Trem- 
blay, M. Civelli and G. Kotliar, Short-range correlation 
induced pseudogap in doped Mott insulators (2005), cond- 
mat/0502565, Phys. Rev. B in press. 

68 D. Senechal, in High Performance Computing Systems 
and Applications and OSCAR symposium, edited by 
D. Senechal (2003). 

69 K. Aryanpour, T. A. Maier and M. Jarrell, Phys. Rev. 
B 71, 037101 (2005). 

70 G. Biroli and G. Kotliar, Phys. Rev. B 65, 155112 (2002). 

71 G. Biroli and G. Kotliar, Phys. Rev. B 71, 037102 

(2005) ; B. Kyung, G. Kotliar, A. -M. S. Tremblay, 
Quantum Monte Carlo Study of Strongly Correlated Elec- 
trons: Cellular Dynamical Mean-Field Theory, (2006) 
cond-mat/0601271. 

72 J. Favand, S. Haas, K. Penc, F. Mila and E. Dagotto, 
Phys. Rev. B 55, R4859 (1997). 

73 C. Dahnken, M. Aichhorn, W. Hanke, E. Arrigoni and 
M. Potthoff, Phys. Rev. B 70, 245110 (2004). 

74 M. Capone, M. Civelli, S. S. Kancharla, C. Castellani 
and G. Kotliar, Phys. Rev. B 69, 195105 (2004). 

75 T. A. Maier, O. Gonzalez, M. Jarrell and T. Schulthess, 
Two quantum cluster approximations (2002), cond- 
mat/0205460. 

76 K. Pozgajcic, Quantitative aspects of the dynamical im- 
purity approach (2004), cond-mat/0407172. 

77 D. Senechal and A.-M. S. Tremblay, Phys. Rev. Lett. 92, 
126401 (2004). 

78 C. Huscroft, M. Jarrell, T. Maier, S. Moukouri and 
A. Tahvildarzadeh, Phys. Rev. Lett. 86, 139 (2001). 

79 C. Grober, R. Eder and W. Hanke, Phys. Rev. B 62, 
4336 (2000). 

80 A. Moreo, S. Haas, A. Sandvik and E. Dagotto, Phys. 
Rev. B 51, 12045 (1995). 

81 S. Moukouri and M. Jarrell, Phys. Rev. Lett. 87, 167010 



(2001). 

82 A.Brooks Harris and Robert V. Lange, Phys. Rev. 157, 
295 (1967). 

83 M.B.J. Meinders, G.A. Sawatzky and G.A. Sawatzky, 
Phys. Rev. B 48, 3916 (1993). 

84 B. Kyung, J. Landry, D. Poulin and A.-M.S. Tremblay, 
Phys. Rev. Lett. 90, 099702 (2003). 

85 F. Ronning et al., Phys. Rev. B 67, 165101 (2003). 

86 N. Armitage et al., Phys. Rev. Lett. 88, 257001 (2002). 

87 B. Kyung, V. Hankevych, A.-M. Dare and A.-M. S. 
Tremblay, Phys. Rev. Lett. 93, 147004 (2004). 

88 Q. Yuan, F. Yuan and C. S. Ting, Phys. Rev. B 72, 
054504 (2005). 

89 H. Kusunose and T. M. Rice, Phys. Rev. Lett. 91, 186407 

(2003) . 

90 C. Kusko, R. S. Markiewicz, M. Lindroos and A. Bansil, 
Phys. Rev. B 66, 140513 (2002). 

91 K. Borejsza and N. Dupuis, Phys. Rev. B 69, 085119 

(2004) . 

92 J. Carbotte, E. Schachinger and D. Basov, Nature 401, 
354 (1999). 

93 D. Scalapino, Physics Reports 250, 329 (1995). 

94 B. Tobijaszewska and R. Micnas, Physica Status Solidi 
B 242, 468 (2005). 

95 A. Moreo and D. J. Scalapino, Phys. Rev. B 43, 8211 
(1991). 

96 R. Scalettar and R. Singh, Phys. Rev. Lett. 67, 370 
(1991). 

97 A. Paramekanti, M. Randeria and N. Trivedi, Phys. 
Rev. Lett. 87, 217002 (2001). 

98 S. Sorella, G. B. Martins, F. Becca, C. Gazza, L. Capri- 
otti, A. Parola and E. Dagotto, Phys. Rev. Lett. 88, 
117002 (2002). 

99 W. Hanke, M. Aichhorn, E. Arrigoni and M. Potthoff, 
Correlated band structure and the ground-state phase di- 
agram in high-T c cuprates (2005), cond-mat/0506364. 

100 C. Pao and N. Bickers, Phys. Rev. B 49, 1586 (1994). 

101 J. Deisz, D. Hess and J. Serene, Phys. Rev. Lett. 76, 
1312 (1996). 

102 F. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

103 O. Andersen, A. Liechtenstein, O. Jepsen and 
F. Paulsen, Journal of the Physics and Chemistry of 
Solids 56, 1573 (1995). 

104 A. Macridin, T. A. Maier, M. S. Jarrell and G. Sawatzky, 
Phys. Rev. B 71, 134527 (2005). 

105 R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, 
C. D. Frost, T. E. Mason, S.-W. Cheong and Z. Fisk, 
Phys. Rev. Lett. 86, 5377 (2001). 

106 J.-Y. Delannoy, Ph.D. thesis, Ecole normale superieure de 
Lyon (2005). 

107 J.-Y. P. Delannoy, M. J. P. Gingras, P. C. W. Holdsworth 
and A.-M. S. Tremblay, Phys. Rev. B 72, 115114 (2005). 

108 L. Raymond, G. Albinet and A.-M. Tremblay, cond- 
mat/0510811. 

109 A. M. Toader, J. P. Goff, M. Roger, N. Shannon, J. R. 
Stewart and M. Enderle, Phys. Rev. Lett. 94, 197202 

(2005) . 

110 A. Damascelli, Z. Hussain and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

111 A. Damascelli, Probing the low-energy electronic struc- 
ture of complex systems by ARPES (2003), cond- 
mat/0307085. 

112 A. Kordyuk and S. V. Borisenko, ARPES on HTSC: 



32 



simplicity vs. complexity (2005), cond-mat/0510218. 

113 p p reusS] \\/ Hanke and W. von der Linden, Phys. 
Rev. Lett. 75, 1344 (1995). 

114 S. Pairault, D. Senechal and A.-M. Tremblay, Eur. Phys. 
J. B (France) 16, 85 (2000). 

115 S. Pairault, D. Senechal and A.-M. S. Tremblay, Phys. 
Rev. Lett. 80, 5389 (1998). 

116 R. Gooding, K. Vos and P. Leung, Phys. Rev. B 50, 
12866 (1994). 

117 K. Shen, F. Ronning, D. Lu, F. Baumberger, N. In- 
gle, W. Lee, W. Meevasana, Y. Kohsaka, M. Azuma, 
M. Takano et al., Science 307, 901 (2005). 

118 T. Yoshida, X. J. Zhou, K. Tanaka, W. L. Yang, Z. Hus- 
sain, Z.-X. Shen, A. Fujimori, S. Komiya, Y. Ando, 

H. Eisaki et al., cond-mat/0510608. 

119 X. Zhou, T. Yoshida, A. Lanzara, P. V. Bogdanov, 
S. A. Kellar, K. M. Shen, W. L. Yang, F. Ronning, 
T. Sasagawa, T. Kakeshita et al., Nature 423, 398 
(2003). 

12(1 T. A. Maier, M. Jarrell, A. Macridin and C. Slezak, 
Phys. Rev. Lett. 92, 027005 (2004); T. Koretsune and 
M. Ogata, J. Phys. Soc. Japan, 74, 1390 (2005). 

121 A. D. Gromko, A. V. Fedorov, Y.-D. Chuang, J. D. Ko- 
ralek, Y. Aiura, Y. Yamaguchi, K. Oka, Y. Ando and 
D. S. Dessau, Phys. Rev. B 68, 174520 (2003). 

122 N. P. Armitage, D. H. Lu, C. Kim, A. Damascelli, K. M. 
Shen, F. Ronning, D. L. Feng, P. Bogdanov, X. J. Zhou, 
W. L. Yang et al., Phys. Rev. B 68, 064517 (2003). 

123 V. Hankevych, B. Kyung, A.-M. Dare, D. Senechal and 
A.-M. Tremblay, in Proceedings of SNS2004 ( 2005), J. 
Chem. Phys. Sol., 67, 189 (2006). 

124 P. K. Mang, O. P. Vajk, A. Arvanitaki, J. W. Lynn and 
M. Greven, Phys. Rev. Lett. 93, 027002 (2004). 

125 M. Matsuda, Y. Endoh, K. Yamada, H. Kojima, 

I. Tanaka, R. J. Birgeneau, M. A. Kastner and G. Shi- 
rane, Phys. Rev. B 45, 12548 (1992). 

126 Y. Onose, Y. Taguchi, K. Ishizaka and Y. Tokura, 
Phys. Rev. B. 69, 024504 (2004). 

127 Y. Onose, Y. Taguchi, K. Ishizaka and Y. Tokura, 
Phys. Rev. Lett. 87, 217001 (2001). 

128 H. Matsui, K. Terashima, T. Sato, T. Takahashi, S.-C. 
Wang, H.-B. Yang, H. Ding, T. Uefuji and K. Yamada, 
Phys. Rev. Lett. 94, 047005 (2005). 

129 R. P. S. M. Lobo, N. Bontemps, A. Zimmers, Y. Da- 
gan, R. Greene, P. Fournier, C. Homes and A. Millis, 
unpublished. 

130 A.J. Millis, A. Zimmers, R.P.S.M. Lobo, N. Bontemps, 

C. C. Homes Phys. Rev. B 72 224517 (2005). 

131 T. Hanaguri, C. Lupien, Y. Kohsaka, D. Lee, 
M. Azuma, M. Takano, H. Takagi and J. Davis, Nature 
430, 1001 (2004). 

132 C. Stock, W. J. L. Buyers, R. Liang, D. Peets, Z. Tun, 

D. Bonn, W. N. Hardy and R. J. Birgeneau, Phys. Rev. 
B 69, 014502 (2004). 

133 M. Norman, D. Pines and C. Kallin, Adv. Phys. 54, 
715 (2005). 

134 P. Fournier, P. Mohanty, E. Maiser, S. Darzens, 
T. Venkatesan, C. Lobb, G. Czjzek, R. Webb and 
R. Greene, Phys. Rev. Lett. 81, 4720 (1998). 

135 E. Motoyama, P. Mang, D. Petitgrand, G. Yu, O. Vajk, 
I. Vishik and M. Greven, unpublished. 

136 J. Sonier, K. Poon, G. Luke, P. Kyriakou, R. Miller, 
R. Liang, C. Wiebe, P. Fournier and R. Greene, Phys. 
Rev. Lett. 91, 147002 (2003). 



137 J. Sonier, K. Poon, G. Luke, P. Kyriakou, R. Miller, 
R. Liang, C. Wiebe, P. Fournier and R. Greene, Physica 
C: Superconductivity and its Applications 408-410, 783 

(2004) . 

138 S. Wakimoto, R. J. Birgeneau, M. A. Kastner, Y. S. 
Lee, R. Erwin, P. M. Gehring, S. H. Lee, M. Fujita, 
K. Yamada, Y. Endoh et al., Phys. Rev. B 61, 3699 
(2000). 

139 S. Wakimoto, J. M. Tranquada, T. Ono, K. M. Kojima, 
S. Uchida, S.-H. Lee, P. M. Gehring and R. J. Birgeneau, 
Phys. Rev. B 64, 174505 (2001). 

140 T. A. Maier, M. S. Jarrell and D. J. Scalapino, Phys. Rev. 
Lett. 96, 047005 (2006). 

141 M. Maple, MRS Bulletin 15, 60 (1990). 

142 R. Haslinger and A. Chubukov, Phys. Rev. B 67, 140504 
(2003). 

143 F. Marsiglio, Phys. Rev. B 73, 064507 (2006). 

144 G. Deutscher, A. F. Santander-Syro and N. Bontemps, 
Phys. Rev. B 72, 092504 (2005). 

145 H. Molegraaf, C. Presura, D. van der Marel, P. Kes 
and M. Li, Science 295, 2239 (2002). 

146 A. Santander-Syro, R. Lobo, N. Bontemps, Z. Kon- 
statinovic, Z. Li and H. Raffy, Europhys. Lett. 62, 568 
(2003). 

147 M. Norman, M. Randeria, B. Janko and J. Campuzano, 
Phys. Rev. B 61, 14742 (2000). 

148 B. Kyung, A. Georges and A.-M. Tremblay, Potential- 
energy (BCS) to kinetic- energy (BEC)-driven pairing in 
the attractive Hubbard model (2005), cond-mat/0508645. 

149 M. Sutherland, D. G. Hawthorn, R. W. Hill, F. Ron- 
ning, S. Wakimoto, H. Zhang, C. Proust, E. Boaknin, 
C. Lupien, L. Tailleler et al., Phys. Rev. B 67, 174520 
(2003). 

150 P. Anderson, Present status of the theory of high-T c 
cuprates (2005), cond-mat/05 10053 published in this is- 
sue of Journal of Low Temperature Physics. 

151 N.E. Bickers, D.J. Scalapino and R.T. Scalettar, Int. J. 
Mod. Phys. B 1, 687 (1987). 

152 S. Allen, B. Kyung, D. Bergeron, V. Hankevych and 

A. -M. Tremblay, unpublished. 

153 D. M. Broun, P. J. Turner, W. A. Huttema, S. Ozcan, 

B. Morgan, L. R., W. N. Hardy and D. A. Bonn, In-plane 
superfluid density of highly underdoped YBa2 Cuz 0§+ x 

(2005) , cond- mat/0509223. 

154 M. Hettler, M. Mukherjee, M. Jarrell and H. Krishna- 
murthy, Phys. Rev. B 61, 12739 (2000). 

lo ° R. Blankenbecler, D. J. Scalapino and R. L. Sugar, Phys. 
Rev. D 24, 2278 (1981). 

156 M. Potthoff, Eur. Phys. J. B (France) 32, 429 (2003). 

157 This additional constraint leads to a degree of consistency 
between one- and two-particle quantities that is absent 
from the standard Hartree-Fock factorization. 

158 rp^g cons ^ an ^ self-energy is not equal to Un- a as it would 
in the trivial Hartree-Fock factorization (see previous 
note). 

159 Th. Maier, M. Potthoff and D. Senechal, unpublished. 

160 The latter result is disputed by a recent DCA calculation 
by Alexandru Macridin, Mark Jarrell, Thomas Maier, P. 
R. C. Kent, cond-mat/0509166. 

161 Kyle Shen, Private communication. 

162 M. Greven, private communication. 

163 I. Herbut, private communication 



